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THE Swift BAT PERSPECTIVE ON NON-THERMAL EMISSION IN HIFLUGCS GALAXY CLUSTERS 
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ABSTRACT 

The search for diffuse non-thermal, inverse Compton (IC) emission from galaxy clusters at hard X- 
ray energies has been underway for many years, with most detections being either of low significance 
or controversial. In this work, we investigate 14-195 keV spectra from the Swift BAT all sky survey for 
evidence of non-thermal excess emission above the exponentially decreasing tail of thermal emission 
in the flux-limited HIFLUGCS sample. To account for the thermal flux contribution at BAT energies, 
XMM-Newton EPIC spectra are extracted from coincident spatial regions so that both the thermal and 
non-thermal spectral components can be determined simultaneously. We find marginally significant 
IC components in 6 clusters, though after closer inspection and consideration of systematic errors we 
are unable to claim a clear detection in any of them. The spectra of all clusters are also summed to 
enhance a cumulative non-thermal signal not quite detectable in individual clusters. After constructing 
a model based on single temperature fits to the XMM-Newton data alone, we see no significant excess 
emission above that predicted by the thermal model determined at soft energies. This result also 
holds for the summed spectra of various subgroups, except for the subsample of clusters with diffuse 
radio emission. For clusters hosting a diffuse radio halo, relic, or mini-halo, non-thermal emission is 
initially detected at the ~ 5ct confidence level, but modeling and systematic uncertainties ultimately 
degrade this significance. This marginal detection is driven by the mini-halo subgroup, suggesting low 
average magnetic field strengths (B ~ 0.1/iG) in the cores of these clusters. 

Subject headings: galaxies: clusters: general — intergalactic medium — magnetic fields — radiation 
mechanisms: non-thermal — X-rays: galaxies: clusters 



1. INTRODUCTION 

A number of observations, mainly at radio frequen- 
cies, have established that relativistic particles and mag- 
netic fields are part o f the intracluster mediu m (ICM) of 
galaxy clusters (e.g.. iGovoni fc Ferettil 120041 ). The large 
(~Mpc) scale, diffuse structures known as radio halos 
and relics are produced by relativistic electrons spiral- 
ing around ~/xG magnetic fields. Because halos and 
relics are not detected in every cluster, but are only 
found in clu s ters with ongoing ma jor merger activity 
(jBuotd l200l ISchuecker et all 12001 . mergers probably 
temporaril y reaccelerate underlying relativistic p opula- 
tions (e.g.. [SarazinlH99l iBrunetti fe Blasil I2005T ). It is 
important to fully characterize the non-thermal phase 
if the dynamics and general state of the ICM is to be 
understood; the proportion of energy tied up in these 
relativistic components, if significant, may bias inferred 
mass estimates necessary to use clusters as cosmologi- 
cal probes (e g JMantz et al.l l2008 HVikhlinin et~aTl l2009l; 
iVanderlinde et al.1 120101 ) . Unfortunately, synchrotron 
emission alone cannot separately determine particle and 
magnetic field energy densities, and so the total energy in 
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the non-thermal phase remains relatively unconstrained. 
However, the electron population can be independently 
observed through inverse Compton (IC) emission due 
to scattering of the ubiquitous Cosmic Microwave Back- 
ground (CMB) photons, which are up-scattered to X-ray 
energies and may be obs ervable if the e lectron popula- 
tion is sufficiently large (jRephaelil 119790 . Detections of 
IC emission, therefore, have the potential to determine 
whether the non-thermal phase is energetically negligible 
or, particularly if the average magnetic field is large, it 
is sizable enough to affect the dynamics and structure of 
the thermal gas. 

Thermal emission clearly dominates at ^keV energies, 
so searches for excess emission due to an IC spectral 
component arc more easily undertaken at very soft or 
hard (> 10 keV) energies. The latter range is partic- 
ularly promising, given the exponential decline in the 
thermal spectrum and the lack of Galactic and solar wind 
charge exchan ge foregrounds that can hamper searches at 
soft energies (iKoutroumpa et aT1l2009l ; ITakei et alj|2007t 
iBonamente et al.ll2009f). In parti cular, the Swift BAT 
all sky survey (jTueller et al.ll2010D provides a deep map 
of hard energy (14-195 keV) emission from which non- 
thermal excesses can be identified. Its uniform coverage 
and impressive sensitivity makes it the most complete 
dataset from whi ch to study the bri ghtest objects in a 
given class (e.g., iWinter et al.l 120091) . Whereas previ- 
ous searches have concentrated on long pointed obser- 
vations of individual clusters, this survey allows a larger, 
more uniform sample to b e searched, as similarly done by 
lAiello et all ((20091 l2010l ) for detected BAT clusters. To 
take full advantage of this capability, we have chosen the 



2 



Wik et al. 



flux-l imited HIFLUGCS sample ([Reiprich fc Bohringcr 
2002;. which contains the brightest clusters in the sky 
outside the Galactic plane. The selection of the bright- 
est clusters may provide the greatest opportunity to de- 
tect IC emission, as in most models the nearest and most 
luminous clusters are expected to have the strongest IC 
signal. Also, because these clusters are bright and con- 
tained within a well-defined survey, there already exist 
good observations at lower X-ray energies, which can be 
used to strongly constrain the thermal properties of the 
ICM - an important prerequisite for the robust detection 
of an IC excess. Finally, the fact that HIFLUGCS is a 
complete flux-limited survey allows one to discuss the 
statistical properties of their hard excesses by stacking 
the individual cluster observations. 

Because they are nearby and bright, many of the clus- 
ters in HIFLUGCS have been targets o f IC searches with 
other telescopes, including A3667 (iFinoguenov et all 
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12005ft. A1367 dHenriksen fc Mushot zkvl I2001D. A2199 
(iKenroner fc Sarazinl 12000ft . and A2163 (|Rephaeli et al.l 
12006ft . Most often clusters are targeted because they 
host a radio halo or relic, as the IC flux then leads to a 
direct measure of the average magnetic field strength. 
A large fraction of HIFLUGCS clusters were also 
included in an analysis of all long expo sure Beppo-SAX 
observations (Ncvalai nen et al.l 12004ft . which found 
marginal evidence for non-thermal excesses in individual 
clusters but a substantial excess in a stacked spectrum. 
In general, an IC component distinct from thermal 
emission in the hard band has been difficult to clearly 
identify, with perhaps the only counter example being 
an exc eptionally deep o bservation of the Ophiuchus 
cluster (|Eckert et al.ll2008ft . The cluster most thoroughly 
searched for non-thermal emission, also in HIFLUGCS, 
is the Coma cluster. Controversial (TRossetti fc Molendi 
detections with RXTE (iRephaeli fc Grubei 



thermal emission; this is particularly true given that the 
Swift BAT survey spectra are coarsely binned (8 channels 
spanning 14 keV < E < 195 keV). Also, any non-thermal 
component in the BAT spectra must be consistent with 
the spectra at softer energies. XMM-Newton is the ideal 
observatory to provide such complementary spectra. For 
one, its large field of view (FOV) allows a higher frac- 
tion of the total emission, which can be quite extended 
given the low redshifts of the sample, to be detected in a 
single pointing. Additionally, the EPIC instruments are 
sensitive to 5-10 keV photons, which make them more 
useful for constraining the highest temperature gas, and 
the telescopes have good spatial resolution so that point 
sources can be excluded from the spectra. Last, but of 
no less importance, XMM-Newton has observed all but 
one (Abell 2244) of the clusters in HIFLUGCS. Unfortu- 
nately, another 4 cluster observations (Abell 401, Abell 
478, Abell 1736, and Abell 2163) are heavily contami- 
nated by background flares and conseq uently unusable 
(for more details, see lZhang et al.ll2011ft . However, the 
data for the remaining 59 clusters are of sufficient qual- 
ity to help constrain potential non-thermal signals in the 
BAT energy bands. 

We extract XMM-Newton spectra for each cluster from 
the largest circular region that either covers the FOV 
or extends to the point where cosmic X-ray background 
(CXB) emission begin s to domina t e, by summing the 
annular spectra from iZhang et al.1 (|2009ft . To ensure 
near Gaussian statistics for \ 2 fitting, adjacent chan- 
nels are grouped until each new bin contains at least 
30 counts. The centers and radii of the circular re- 
gions, along with each pointing's observation ID, are 
listed in Table [TJ Source spectra are extracted in con- 
centric annuli within the region; corresponding parti- 
cle background spectra are derived from CLOSED mode 
calibration data, which are renormalizcd based on 3-10 
keV events out of the FOV and outside of a 15'. 4 ra- 
di us from the detecto r center (for details see Section 2.4 
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have recently be en challenged with comparable Suzaku 
(|Wik et al.ll2009h observatio ns and a detailed analysis of 
the Swift BAT survey data (|Wik et al.ll2"oTTI ). 

To perform the deepest hard X-ray survey of non- 
thermal emission in clusters to date, we jointly fit high 
quality XMM-Newton EPIC and Swift BAT spectra, ex- 
tracted from identical regions and cross-calibrated to 
make their absolute spectral responses as consistent as 
possible. We describe the data and its calibration in 
Section [5] In Section [31 the thermal and non-thermal 
character of the spectra are separately analyzed, and in 
Section [4] they are jointly fit for each individual clus- 
ter. We also search for a statistical hard excess in sets 
of stacked spectra for the entire sample and for several 
subsamples in Section [5l Lastly, the implications of our 
results are discussed in Section [5] We assume a flat cos- 
mology with n,M = 0.23 and Hq = 70 km s _1 Mpc -1 . 
Unless otherwise stated, all uncertainties are given at the 
90% confidence level. 

2. OBSERVATIONS AND DATA PREPARATION 

2.1. XMM-Newton EPIC Spectra 

For the lower energy BAT bands, it is very useful to 
have X-ray spectra at lower energies to constrain the 



is described in IZhang et al.l (|2009f ) . As an additional 
step, we readjust the normalization of the particle back- 
ground spectra by hand to ensure the 7-12 keV contin- 
uum of the cluster spectra have a more physical shape. 
We define "more physical" as the background normal- 
ization that minimizes the % 2 statistic for a single tem- 
perature (IT) (using the APEC plasma emission modeQ) 
individually fit to the EPIC-pn (2 < E < 12 keV) and 
MOS1 and MOS2 (2 < E < 10 keV) spectra. The new 
best-fit temperatures, after these initial renormalizations 
of the background, are compared to each other and to 
previ ous measurements (primarily iReiprich fc Bohringerl 
2002). While this method may bias the background level, 
especially if a single temperature model is a poor descrip- 
tion of a given spectrum, repeating this procedure with 
two temperature (2T) and single temperature plus power 
law (T+NT) models yield comparable or inferior results, 
usually favoring obvious under-subtractions of the back- 
ground that produce systematic patterns in the residu- 
als. We favor normalizations that leave the background 
slightly under-subtracted, in order to avoid removing a 
real non-thermal signature. For the most part, the over- 
all spectrum is only mildly affected since much of the 

7 http:/ /cxc. harvard.edu/atomdb/sources_apec. html 
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emission is at lower energies where the background is a 
smaller fraction of the total. One consequence is that 
instrumental lines, which are typically between 7.5 and 
9.5 keV and mainly are a problem in the EPIC-pn spec- 
tra and which can vary in intensity relative to the back- 
ground continuum, can be under- or over-subtracted. No 
resolved ICM lines exist in this range, so we simply ignore 
th is energy range w hen poor line subtractions occur, as 
in iWik et al.l (|2009l ). Based on the change in x 2 as the 
background normalization is varied, a typical 90% level 
uncertainty in the normalization is ~ 3%. 

We choose to model, instead of subtract, one further 
backgrou nd componen t : the CXB due to extragalactic 
sources. iLumb et al.l (2002), using XMM-Newton sky 
fields, find that this component of the CXB is well fit 
by a power law with photon index of 1.42 in the hard 
band (2-10 keV). Their results are in good agreement 
with other work in this b and (e.g.. iMorctti et al.l 120031 : 
iDe Luca k. Molendl I2004D . We adopt their normaliza- 
tion at 1 keV of 8.44 photons cm -2 s _1 keV -1 sr _1 , 
which is scaled to match the extraction area for each clus- 
ter. The impact of cosmic variance, or the field-to-field 
variation in CXB flux resulting from large scale struc- 
ture and source population selection, is not included as 
a systematic uncertainty in the following analysis due to 
its small effect. While cosmic variance increases with de- 
creasing solid angle, the high sensitivity of XMM-Newton 
allows most of the sources responsible for a higher vari- 
ance to be removed, so for one of our typical regions the 
90% uncertainty i s only -10% of the CXB flux. Note 
that ILumb et al.1 ([2002) remove detected point sources 
as is done here, so their spectrum can be directly applied 
as is. The Galactic component of the CXB is also not 
considered, as it only contributes below 1 keV, and we 
restrict our fits to the 2-12 keV range. 

2.2. Swift BAT 58-month Survey Spectra 

The Swift mission and the prope r ties o f the survey 
are des cribed in detail in IWik et al.l (|2011l Section 2.2) 
and in Tucllcr et al. ([2O1O0- Similarly, we refer to that 
Section and the appendices for details on the extraction 
and calibration of sources from survey image data. To 
briefly summarize, the flux calibration is tied to the Crab 
spectrum, which we define to have the same spectrum 
as that observed by XMM-Newton for E > 2 keV, ex- 
trapolated to BAT energies via an adopted model based 
on Suzaku observations. In this way, both the cross- 
normalization and spectral shape of the XMM-Newton 
and Swift spectra will match, and continuous models can 
be jointly fit to them simultaneously. Unfortunately, 
independent measurements by various instruments, in- 
cluding the Swift BAT, have recently demonstrated that 
the hard X-ray spect rum of the Crab is in fact variable 
on yearly timescales (jWilson- Hodge et"al~ll2011[ ) . At 14- 
50 keV energies, where the only appreciable amount of 
flux is detected from clusters, the variation spans about 
10% over the last 5 years, with a co nsistent decline only 
over the last 2 years (see Fig. 5 of iWilson-Hodge et al.l 
2011). Our flux calibration of BAT sources depends on 
an adopted model for the Crab spectrum, which is taken 
from Suzaku XIS and HXD-PIN observations that took 
place in August 2005; this time occurs during one of 
the higher flux periods. Since the BAT survey spectrum 
of the Crab spans the following 5 years of observations 



TABLE 1 

XMM-Newton Observations, Regions, and Cluster Classes 



a(J2000) <5(J2000) Radius 



Nciinc 


OKcTFl 
\J DS1U 


\ ae S) 


(deg) 


(arcmin) 




A0085 


0065140101 


10.45957 


-9.30303 


11.6667 


SCC 


A0119 


0505211001 


14.07130 


-1.25327 


9.3333 


NCC 


A0133 


0144310101 


15.67971 


-21.87968 


6.7000 


SCC 


NGC507 


0080540101 


20.91068 


33.25063 


9.4667 


SCC 


A0262 


0109980101 


28.19002 


36.15114 


13.5333 


SCC 


A0400 


0404010101 


44.42226 


6.02696 


12.6667 


NCC 


A0399 


0112260101 


44.46513 


13.04713 


10.4000 


NCC,R 


A3112 


0105660101 


49.49456 


-44.23562 


6.9667 


SCC 


Fornax 


0400620101 


54.61989 


-35.45122 


10.6333 


SCC 


2A0335 


0109870101 


54.66787 


9.96803 


8.9667 


SCC.R 


IIIZw54 


0505230401 


55.32801 


15.40390 


6.8667 


WCC 


A3158 


0300211301 


55.72316 


-53.63099 


9.1333 


NCC 


NGC1550 


0152150101 


64.90839 


2.40929 


11.6667 


SCC 


EXO0422 


0300210401 


66.46339 


-8.56118 


7.1333 


SCC 


A3266 


0105260901 


67.81198 


-61.44835 


12.0000 


WCC 


A0496 


0135120201 


68.40753 


-13.26069 


10.1667 


SCC 


A3376 


0151900101 


90.54203 


-39.95994 


6.0000 


NCC,R 


A3391 


0505210401 


96.60081 


-53.69002 


6.7333 


NCC 


A3395s 


0400010301 


96.69188 


-54.54530 


4.2000 


NCC 


R1504 


0401040101 


106.37174 


-12.93125 


8.1992 


SCC,R 


A0576 


0205070301 


110.35886 


55.75948 


9.0000 


WCC 


A0754 


0136740101 


137.32574 


-9.68781 


10.1667 


NCC,R 


HydraA 


0109980301 


139.52491 


-12.09342 


5.0000 


NCC 


A1060 


0206230101 


159.17853 


-27.52841 


8.4667 


WCC 


A1367 


0061740101 


176.18539 


19.73211 


10.0000 


NCC,R 


MKW4 


0093060101 


181.11522 


1.89480 


8.3333 


SCC 


ZwC11215 


0300211401 


184.41928 


3.65818 


6.2333 


NCC 


NGC4636 


0111190701 


190.70940 


2.69179 


9.8333 


SCC 


A3526 


0406200101 


192.21101 


-41.30430 


12.9333 


SCC 


A1644 


0010420201 


194.29469 


-17.40291 


14.7333 


SCC 


A1650 


0093200101 


194.67448 


-1.75920 


5.1667 


WCC 


A1651 


0203020101 


194.84310 


-4.19633 


7.5000 


WCC 


Coma 


0124711401 


194.93888 


27.95150 


14.6667 


NCC,R 


NGC5044 


0037950101 


198.84908 


-16.38664 


11.5000 


SCC 


A3558 


0107260101 


202.00169 


-31.50027 


10.2333 


WCC 


A3562 


0105261801 


203.40201 


-31.67382 


6.1667 


WCC,R 


A3571 


0086950201 


206.86609 


-32.86052 


8.5000 


WCC 


A 1795 


0097820101 


207.21991 


26.59282 


8.0000 


SCC 


A3581 


0205990101 


211.87760 


-27.01320 


11.0667 


SCC 


MKW8 


0300210701 


220.17560 


3.47159 


7.5667 


NCC 


A2029 


0111270201 


227.73326 


5.74264 


6.5000 


SCC.R 


A2052 


0109920101 


229.18501 


7.02012 


7.0667 


SCC 


MKW3S 


0109930101 


230.45945 


7.70323 


9.6667 


SCC 


A2065 


0112240201 


230.62112 


27.72063 


6.6667 


WCC 


A2063 


0550360101 


230.77401 


8.60701 


7.1667 


WCC 


A2142 


0111870301 


239.56451 


27.25178 


6.6667 


WCC,R 


A2147 


0505210601 


240.56789 


15.97177 


11.3333 


NCC 


A2199 


0008030201 


247.15461 


39 5481 1 

'-J \J ■ Ortt^J J. _L 


12 3333 


SCC 


A2204 


0112230301 


248.19604 


5.57554 


6.1333 


SCC,R 


A2256 


0141380201 


255.96829 


78.67197 


8.0000 


NCC,R 


A2255 


0112260801 


258.22709 


64.06428 


8.1667 


NCC,R 


A3667 


0206850101 


303.16966 


-56.84081 


13.0000 


WCC,R 


S1101 


0123900101 


348.49294 


-42.72664 


6.0333 


SCC 


A2589 


0204180101 


350.98652 


16.77595 


5.0000 


WCC 


A2597 


0147330101 


351.33334 


-12.12416 


6.5667 


SCC 


A2634 


0002960101 


354.62099 


27.03107 


11.0000 


WCC 


A2657 


0402190301 


356.23640 


9.19810 


5.6667 


WCC 


A4038 


0204460101 


356.93602 


-28.14506 


12.3333 


WCC 


A4059 


0109950201 


359.25704 


-34.75803 


9.1333 


SCC 


From Hudson et al.1 (|2010T): SCC 


= "strong 


cool core 


cluster," 



WCC = "weak cool core cluster," and NCC = "non-cool core clus- 



ter"; clusters that host a radio halo and/or relic are labeled with 
"R" 

and averages over these fluctuations, the normalization 
of our adopted model is only 2-3% higher than the ac- 
tual flux emitted. The effect this has on our derived 
fluxes is to make them 2-3% higher than they actually 
are; this amount is equivalent to the la error on the 
14-20 keV flux of Coma, which is the highest signal- 
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to-noise flux consid e red h ere by a factor of 2. Also, 
iWilson-Hodge et all (|201lD show that the recent decline 
in flux is more dramatic for higher energy bands. Both 
of these behaviors - the overall decline in flux and the 
steepening of the spectrum - bias our derived BAT fluxes 
high, which could lead to a higher chance of false non- 
thermal detections. However, since we find no convincing 
evidence for non-thermal excesses even given this proba- 
ble effect, and since we allow for a 10% cross-calibration 
uncertainty between the BAT and EPIC spectra, which 
easily encompasses this level of variability, it is clear that 
our choice of flux calibrator does not strongly impact 
the following analysis, except to make our upper limits 
slightly more conservative than they otherwise would be. 

While the standard processing of coded mask imaging 
data is designed to extract the fluxes of point sources, 
it is also possible to extract the flux of a mildly ex- 
tended source, albeit with somewhat g reater uncertainty 
([Renaud et all [20061: iWik et al.l l201lW . The large effec- 
tive PSF (full width at half maximum FWHM ~ 20') 
for point sources in the survey means that even nearby 
clusters of galaxies will appear only slightly extended; the 
FWHM of the Coma cluster - the most extended, reliably 
detected source in the survey - is only 28.'5. Note that 
while 4 clusters (Fornax, NGC 4636, A3526, and A1060) 
have larger ang ular extents than Co ma (based on angular 
i?50o estimates. lEckert et al.l (120111 )1. they are all cooler, 
less massive systems and thus either not detected or only 
marginally detected by the BAT at 14-20 keV. From Fig- 
ure [TJ it is clear that detected clusters (colored circles) 
are typically extended, relative to other sources. The 
horizontal lines mark the standard deviation of best-fit 
FWHM values for the non-cluster sources in each signal- 
to-noise bin; they also represent the approximate error 
on FWHM estimates for the clusters in each bin. Indi- 
vidual clusters are labeled in the 4 lowest energy BAT 
bands when they are detected at a signal-to-noise ra- 
tio greatei_than_5. We follow the procedure outlined 
in lWik et al.l (|2011l ) to extract fluxes for diffuse sources, 
which requires the spatial distribution of the emission 
to be known. Because clusters are comparable in size to 
the effective spatial resolution of the survey, detailed spa- 
tial models are not necessary to extract accurate fluxes. 
We consider generic /3-model surface brightness profiles, 
which well represent the radial profiles at softer energies. 
Taking a representative value for /3 of 0.75, we find that 
all > 3er detected clusters (in a given band) can be well fit 
with core radii r c of either 4', 6', 8', or 10'. Profiles with 
r c < 4' are hard to distinguish from point source profiles, 
so for any cluster emission that is too narrow to be fit 
with the r c = 4' model is treated as a point source. The 
true spatial distribution may differ from these fiducial 
models, but our aim is only to extract accurate fluxes, 
not describe the distribution of hard X-ray emission. For 
Coma, a /3-model fit in the first BAT band (El: 14-20 
keV) yields a total flux 9% lower than that derived from 
a more detailed model of its spatial distributi on derived 
from an XMM-Newton temperature map (see IWik et aD 
1201 lh . which accounts for th e NE-SW non-axisymmetric 
elongation of the emission (jEckert et al.l [2001 . While 
9% is a significant difference, Coma is one of the most 
significantly detected and is the most extended cluster 
in the survey, so this deviation, which amounts to a fac- 
tor of only 1.6 times the 1-a error on the flux, is the 



largest we would expect using this set of extended mod- 
els. Also, note that no energy dependence in FWHM 
values is detected; e.g., Coma shows some variation with 
energy band, but these measurements are all consistent 
within their uncertainties. 

We also investigated the use of diffuse models for all 
the clusters, irrespective of their observed extent, to ac- 
count for the possibility that we are missing low surface 
brightness emission obscured by noise. Since the spa- 
tial distribution of E > 10 keV emission is unknown, 
we as sume /3-model profiles deri ved from ROSAT im- 
ages (jRciprich fc Bohringerl [2002h . For clusters with a 
clearly extended BAT profile, these models reasonably, 
but usually not perfectly, follow the emission; however, 
these profiles cannot be reliably distinguished from those 
at lower energies given that background fluctuations can 
still distort the profile due to the low signal-to-noise ra- 
tios. Spectral fits using these fluxes produce similar re- 
sults to those we present in this work, but because their 
associated errors are larger, these spectra are generally 
less sensitive, so any additional flux captured — which is 
not significant — is also diluted. Therefore, these spectra 
are not considered further. 

For clusters with modeled extended emission, we do 
not want to include the portion of flux that falls out- 
side the XMM-Newton extraction region during joint fits 
of the data, since the complementary softer flux in the 
XMM-Newton band spectra is not present. Therefore, 
only the fraction of the flux that resides within the XMM- 
Newton region is included in the spectra derived here. 
One uncertainty, particularly when emission is detected 
at lower significance, is where the emission is actually 
coming from, given the positional accuracy of the sur- 
vey (a 5a source detected in a given band has a 90% 
error circle of radius 6'). Since the El band-derived 
positions are near the center of the extraction region, 
within their respective error circles, we assume the cen- 
ter of the hard band distribution is coincident with the 
center of the XMM-Newton extraction region except for 
A754, A3266, and A2256. For these detected clusters, 
their BAT positions are somewhat offset from the surface 
brightness peak due to an anisotropic te mperature dis- 
tribu tion produced by mergers (see, e.g.,lHenrv fc Briell 
[19951: IFinoeuenov et al. 2006: S un et al.12 002'). Following 
this procedure, we will not underestimate the coincident 
flux, although overestimates may result that could lead 
to incorrect hard excesses. However, since we are unable 
to significantly detect non-thermal emission individually 
in any of the clusters, this procedure can only cause us 
to be biased in favor of more conservative upper limits. 

3. SEPARATE FITS TO INDIVIDUAL XMM-Newton 
EPIC AND Swift BAT SPECTRA 

Before combining the Swift and XMM-Newton 
datasets, we characterize each telescope's spectra sep- 
arately. The goal is to identify any problems with the 
data or our methodology that might lead to biased re- 
sults when the spectra are fit jointly. 

3.1. Single Temperature Fits to the EPIC Spectra 

The motivation for including XMM-Newton spectra in 
the analysis is to fully characterize the thermal proper- 
ties of the hottest gas in the ICM, which will contribute 
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Fig. 1. — Source extent as a function of approximate signal-to-noise ratio (maximum pixel flux divided by local blank sky fluctuations) 
in the 4 lowest energy BAT bands. Non-cluster sources are shown as small dots (omitted in astro-ph version) and horizontal lines mark the 
standard deviation of best-fit FWHM values for the non-cluster sources in each signal-to-noise bin. Galaxy clusters (colored circles) clearly 
trend above the mean FWHM of ~ 20'; clusters are labeled for S/N values above 5; the only cluster detected at E > 50 keV is Perseus, 
but that emission is primarily due to the AGN in the center of NGC 1275. Many bright galaxy clusters are at least somewhat resolved by 
the BAT, which should be accounted for when extracting fluxes from the BAT survey. However, FWHM estimates for sources below a S/N 
~ 10 can be particularly contaminated by background fluctuations and may not be representative of their true spatial extent; the la error 
on the FWHM estimates is approximately given by the horizontal lines bounding the spread in point source FWHM. 



flux to the BAT energy bands. Similarly, these lower en- 
ergy spectra must be consistent with any indication of a 
non-thermal component in the BAT spectra; for exam- 
ple, a steep power law may best describe the BAT data 
but at lower energies result in a poor description of the 
spectrum. Since our purpose is not to fully character- 
ize the total emission detectable by XMM-Newton, but 
only capture the state of the hottest gas, we ignore all 
events with energies below 2 keV. Cool (< 1 keV) gas 
is completely unimportant at BAT energies, and it will 
not overly bias E > 2 keV data. We therefore initially 
consider EPIC spectra in the 2-12 keV range for the pn 
and 2-10 keV range for the MOS detectors; including 
photons down to 2 keV provides additional leverage dur- 
ing spectral fitting, since most of the detected photons, 
regardless of temperature, are at lower energies. 

However, the lower end of this energy range presents 
two issues. First, bright ~ 1 keV gas can significantly 
contribute to the emission between 2 and 3 keV, which 
certainly exists in some of the cool core clusters in HI- 
FLUGCS. In single temperature fits, the average tem- 
perature will then be biased low to accommodate this 
component, which could lead to thermal emission being 
interpreted as a non-thermal excess. Multi-temperature 
fits would alleviate this problem, but most of the XMM- 
Newton data are not of sufficient quality to strongly con- 



strain more than one temperature component in this en- 
ergy range. Including E < 2 keV data to better con- 
strain multi-temperature fits would also require a more 
complicated analysis that will involve more free param- 
eters and, because the highest signal-to-noise ratios are 
in the ~ 1 keV channels, fits would be driven by this 
data, possibly resulting in biased high temperature com- 
ponents. The second issue relates to the imperfectly cal- 
ibrated gold edge at 2.2 keV, where the response drops 
somewhat abruptly. While on its own this feature does 
not strongly impact spectral fits, because it lies near the 
edge of our energy range where the signal-to-noise ratio is 
largest, secondary model components can be "co-opted" 
into better fitting this edge. For instance, in a spectrum 
truly described by a gas at a single temperature, the addi- 
tion of a second temperature or non-thermal component 
to the fit will cause the second component to "fix" any 
deviations at this edge, typically resulting in a low tem- 
perature or steep photon index that has no real physical 
counterpart. 

In practice, both of these effects can conspire to pro- 
duce the appearance of a more significant non-thermal 
spectral component than is warranted by the rest of the 
data. To counter both issues, we also perform fits to 
data with energies E > 3 keV, which exclude the gold 
edge and any sizable emission from < 1 keV gas. These 
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spectra have lower signal-to-noise due to excluding the 
2-3 keV emission, but the high fluxes of clusters in our 
sample reduce this issue's importance. Single temper- 
ature fits in both the 2-12 keV and 3-12 keV ranges, 
jointly fit to all three EPIC spectra (except for A3526, 
for which the MOS-1 spectrum is ignored, and for A2142 
and A2147, for which the MOS-2 spectra are ignored), 
are given in Table [2j The pn and MOS instrument cross- 
normalization is left as a free parameter, which allows for 
a typ ical (10 ±10)% difference between their calibration 
(e.g.. lSnowdenll2002t ). This cross-normalization factor is 
used and kept fixed during all subsequent joint EPIC- 
BAT fits. The change in the best-fit temperature from 
the E > 2 keV to E > 3 keV fits is only ~ 0.3 keV on 
average, indicating that the temperature is generally ro- 
bust to the choice of the energy range, but that higher 
energy photons come preferentially from higher tempera- 
ture gas, assuming the true temperature structure is not 
isothermal but contains a continuous spectrum with gas 
at many t emperatures due to substructure and/or radi al 
gradients (|Cavagnolo et al.ll2008t fenowden et al.ll2008t) . 

3.2. Non-thermal Fits to the BAT Spectra 

Our goal is to detect a non-thermal spectral compo- 
nent at hard energies, but because the statistical weight 
of the BAT channels is so much less than the EPIC chan- 
nels (lower S /N and fewer of them, at least by an order of 
magnitude), we have to be careful not to let the XMM- 
Newton data unfairly drive the spectral fits. To assess 
the sensitivity of our BAT spectra, we extract 10,000 
blank sky spectra from uniformly distributed, random 
positions at least 40' from any known sources and greater 
than 20° from the Galactic plane, to mimic the selection 
function in HIFLUGCS. We then fit these spectra with 
a fiducial power law model of photon index Y fixed at a 
value of 2, roughly the appropriate slope for IC emission 
inferred from radio halos, relics, and mini-halos. While 
the spectral index determined from the radio is typically 
steeper than this (2.2-2.4), the electrons producing the 
radio emission at v > 100 MHz have higher energies than 
those producing IC at E < 50 keV for B < 0.5/iG, 
so a simple extrapolation may not be appropriate. A 
clear flattening of the radio spectrum a t low frequencies 
is ap parent in so me cases, e.g.. Coma (IThierbach et al.l 
120031 ) and A3562 ([Giacintucci et al.l l2005D. although this 
is not u niversally found as in A2256 (|Brentiensl[2008l ) and 
A2255 (jPizzo fc de Bruvnll2009f ) . Since the BAT data are 
not particularly sensitive to the precise value of the in- 
dex, we choose a natter slope to avoid poorly fitting the 
data at ~ keV energies where the power law distribution 
of relativistic electrons is most likely to turn over in a 
steady state-like injection model (e.g.. lSarazinlll999l ). 

The distribution of best-fit normalizations from these 
power law fits are presented in the narrow histogram in 
Figure [3J They are well fit by a symmetric Gaussian 
(dashed smooth line) and indicate a la sensitivity thresh- 
old of - 2 x 10" 12 erg cm" 2 s" 1 (20-80 keV). Similarly, 
the formal 3<r detection level is 5.8 x 10" 12 erg cm" 2 s" 1 . 
In principle, the BAT survey is sensitive enough to con- 
firm or reject previous detections of hard excesses with 
fluxe s ~ 10~ n erg cm ~ 2 s" 1 (e.g.. iRephaeli fc Gruberl 
[20021: iMoIendi et aij|200l iFusco-Femiano et al.ll2004D . 

Now we wish to compare our cluster spectra with this 
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Fig. 2. — The distribution of power law normalizations (with a 
fixed photon index T = 2) fit to 10,000 blank sky spectra extracted 
from the BAT survey (narrowly-binned histogram). The best-fit 
Gaussian distribution is overlaid as the smooth, dashed line (red). 
Similar best-fit normalizations are shown for the 59 HIFLUGCS 
clusters (see text for details), with individual normalizations rep- 
resented as vertical lines (blue). The cluster histogram (wide bins) 
has been scaled up to show its agreement with the blank sky spec- 
tra. In general, the cluster BAT spectra lack any clear evidence 
for a non-thermal component, except in a few cases comprising the 
positive tail of the blue histogram. 

distribution, but first we have to account for any thermal 
emission in the lower energy bands. The single temper- 
ature models derived with XMM- Newton (2-12 keV) are 
included as a second component along with the power law 
model, with only its normalization left as a free param- 
eter. The resulting non-thermal normalizations are also 
given in Figure [2] as both the wider histogram (scaled 
up) and as the vertical lines (showing individual val- 
ues). While the majority of cluster non-thermal com- 
ponents are consistent with the blank sky fits, there is 
a tail at positive normalizations possibly indicative of a 
non-thermal excess. However, the thermal contribution 
is not well determined in this method and may be un- 
derestimated. Intriguingly, the three clusters with the 
most significant non-thermal component (A2029, A1367, 
and A1651) have positive fluxes, although marginally de- 
tected, in all 8 BAT bands; this rarely occurs for the 
blank sky spectra. We discuss these clusters in more de- 
tail later. The main result from this analysis is that the 
BAT cluster spectra have probably not reached a sen- 
sitivity level sufficient to detect hard, non-thermal ex- 
cesses, if they exist, in the brightest clusters. 

4. JOINT FITS TO THE EPIC-BAT SPECTRA 

BAT fluxes are calibrated to match both the normal- 
ization and the spectral shape of sources a s detected 
by th e XMM-Newton EPIC-pn instrument ([Wik et al.1 
12011 . and they are extracted from regions identical to 
the XMM-Newton extraction regions. As such, contin- 
uous spectral models can be used over the full 2-195 
keV energy range to simultaneously fit both the XMM- 
Newton and Swift spectra. However, in individual cases 
the cross-normalization factor, fcNi may stray from a 
value of 1 as it does between the pn and MOS instru- 
ments (see Section |3~TT) . We therefore adopt, along with a 
3% uncertainty in the XMM-Newton background normal- 
izations, a conservative 10% systematic uncertainty for 
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Fig. 3. — A comparison of best-fit temperature values in IT fits 
to only the XMM-Newton spectra (x-axis) and to the EPIC and 
BAT spectra simultaneously (y-axis). Solid lines indicate the 90% 
error interval for E > 2 keV fits, dashed (red) lines for E > 3 keV, 
and the dashed diagonal line represents equality between the two 
temperature determinations. Jointly fitting both datasets yields 
consistent temperatures to those derived only in the XMM-Newton 
band. Fitting over a slightly higher energy range {E > 3 keV), 
while increasing the average temperature by ~ 0.3 keV in the EPIC 
bandpass (see text), does not increase the joint fit temperatures as 
much; note how the dashed points fall slightly below equality for 
moderately hot clusters. 

fcN- Because no compelling evidence for non-thermal 
emission is found in the nominally calibrated spectra 
(see analysis below), we only consider these uncertain- 
ties when deriving 90% confidence interval upper limits. 

4.1. General Properties from the Joint Analysis 

For each cluster, 3 simple spectral models are employed 
to describe the emission covering 2 orders of magnitude 
in energy: a single temperature thermal model (IT), a 
two temperature model (2T), and a thermal plus non- 
thermal model (T+IC). Due to the limited sensitivity 
of the Swift data, more complicated models cannot be 
constrained; for example, the separate temperature com- 
ponents in the 2T model are generally poorly constrained 
in our analysis. Above 50 keV, the APEC emission model 
is replaced with MeKa because APEC is not defined above 
50 keV in the implementation of XSpec used here (Ver- 
sion 12.6.0k). Note that the MeKaL emission model could 
also be used continuously across this energy range, if 
the look-up table switch is turned off. For the thermal 
component, the temperature, abundance, redshift, and 
normalization are all varied. The individual abundances 
and redshifts in the 2T model are tied together. The 
non-thermal photon index is initially fixed at T = 2, 
typical of radio halos, and the normalization is allowed 
to vary; when the photon index is fit for, it is always 
fixed to the best-fit value before errors for other param- 
eters are derived. In general, the photon index is poorly 
constrained, allowing for a wide range of normalizations, 
which are then less straightforward to evaluate. The pur- 
pose of fitting for the photon index is to make sure that 
we are not biased against detectable IC components with 
indices that differ from the fiducial value. 

Because of complications arising at energies between 2 
and 3 keV (see Section l3Tj) . we perform these fits for both 
the 2-195 keV (Table EJ and the 3-195 keV (Table [BJ 




kT (keV) [1T Model (3-195 keV)] 

Fig. 4. — Upper limits and measurements of the non-thermal 
spectral component in the 3-195 keV joint fits as a function of 
cluster temperature. Limits and error bars indicate the 90% con- 
fidence interval without considering the impact of systematic un- 
certainties. In general, an excess attributable to IC emission is 
not observed, and the few detections, discussed individually in the 
text, have marginal statistical significance. 

spectral ranges. The E > 2 keV fits, at first glance, sug- 
gest that there may be evidence for a non-thermal com- 
ponent in a majority of HIFLUGCS clusters. Many of the 
clusters with some evidence, at least at the 90% level, of 
a non-thermal excess are, unexpectedly, low temperature 
clusters without significant detections at BAT energies. 
In these cases, the non-thermal component is serving to 
"adjust" a problem at lower energies - due to either in- 
completely modeled low temperature components, an im- 
perfectly calibrated response at the gold edge, or both. 
The significance of these instances will disappear from 
fits within a slightly higher energy range, while real non- 
thermal emission will become a higher proportion of the 
total flux and so this component should not greatly di- 
minish in significance. A drastic reduction in the number 
of marginally detected non-thermal excesses is seen when 
comparing Tables [5] and |U only 6 clusters are detected 
to have such emission at the 90% confidence level (sta- 
tistical). These clusters will be discussed individually in 
Section 1431 

While the 3-12 keV band avoids some possible sys- 
tematic uncertainties with the XMM-Newton response 
and complications from cooler gas, the narrower range 
may reduce our ability to strongly constrain multi- 
temperature components in the spectra. One concern 
is that a weak non-thermal emission component might 
be indistinguishable from a purely thermal model with 
a slightly elevated temperature. Note, however, that the 
3-12 keV band temperatures in Section I3TT1 are typically 
only ~ 0.3 keV higher than the 2-12 keV temperatures. 
Therefore, the IT model temperatures should agree for 
the joint fits over both energy ranges, which is found 
to be the case in Figure |3] Temperatures derived from 
joint fits are consistent with those found using only the 
XMM-Newton spectra, for both energy ranges. For the 
most part, temperatures from the joint fit 3-195 keV fits 
are in good agreement with or slightly lower than the 
3-12 keV temperatures. The contribution of the BAT 
data in this case is to somewhat lower the best-fit tem- 
perature, contrary to the expectation if a detectable non- 
thermal excess were present. The 3-195 keV non-thermal 
flux limits and possible detections (90%, statistical) are 
shown in Figure [4] 
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Fig. 5.— Abell 2029: The T+IC model simultaneous fit to the 
EPIC (E < 12 keV) and BAT (E > 14 keV, black data points and 
residuals) spectra. The EPIC-pn spectrum and residuals are in 
green, and the MOS 1 and 2 spectra/residuals are in dark and light 
blue, respectively. The like-colored lines below these spectra show 
the CXB model contribution. The total model fit and thermal 
contribution is represented by the black histogram, and the red 
lines represent the non-thermal (r = 2) spectral component. 
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FlG. 6.— Abell 1651: The T+IC model simultaneous fit to the 
EPIC (E < 12 keV) and BAT (E > 14 keV) spectra. The notation 
is identical to Figure 

4.2. Individual Cases 

Six clusters have a formal detection of non-thermal 
emission in the 3-195 keV band. Two of these 6 clusters 
are also in the top 3 of candidates for emission based on 
their BAT-only fits: A1651 and A2142. The other cluster 
in this top 3 - with the largest non-thermal normaliza- 
tion of all the clusters - is A2029, so we will include this 
cluster with the 6 "detected" clusters as worth some brief 
discussion. The clusters are listed in order of decreasing 
non-thermal flux. 

A2029 (Fig.©: This hot (~ 8 keV), coo l core cluster 
has b een studied in detail with Chandra (Clark e et al.1 
2004) , who explore the interaction between cool gas and 
the radio AGN in the cluster center. The cluster is elon- 
gated but relatively regular; no evidence exists for major 
merger activity; however, a minor merger may be pro- 
ducing the spiral surface brightness enhancement in the 
center. Also, no evidence for an X-ray counterpart of the 
AGN is visible in the Chandra data. In addition to the ra- 
dio jets, the cor e of the cluster is als o host to an extended 
radio minihalo (|Murgia et al.|[2009l) . As with radio halos 
and relics, IC emission may be detec table from th e mini- 
halo if the magnetic field is small; I Taylor et al.l (|1994[ ) 
measured a lower limit of B > 0.11-0.19 /xG with Fara- 
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Fig. 7.— Abell 2142: The T+IC model simultaneous fit to the 
EPIC (E < 12 keV) and BAT (E > 14 keV) spectra. The notation 
is identical to Figure [5] 
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Fig. 8.— Abell 3112: The T+IC model simultaneous fit to the 
EPIC [E < 12 keV) and BAT (E > 14 keV) spectra. The notation 
is identical to Figure [5] 




Energy (keV) 

FlG. 9. — Abell 1367: The T+IC model simultaneous fit to the 
EPIC (E < 12 keV) and BAT (E > 14 keV) spectra. The notation 
is identical to Figure [5] 

day RM observations of the jet. The implied magnetic 
field strength, if we take as the IC flux that found with 
the 2-195 keV fit, is B ~ 0.08 /iG, roughly consistent 
with their field strength. 

But have we really detected IC from the cluster core? 
The significance of the non-thermal component com- 
pletely disappears in the 3-195 keV fit; all three model 
combinations match the data equally well. Also, the 2T 
model formally provides a better fit to the 2-195 keV 
spectrum where the non-thermal component is detected. 
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The second temperature component, ~ 0.3 keV, is con- 
sistent with a low tempe r ature component of 0.11 keV 
observed by iClarke et al.l {2004). Given these results, it 
is more likely that the non-thermal component is try- 
ing to mimic the low kT cool core component in the 
2-3 keV range, since its significance disappears if this 
energy range is ignored. However, it is worth noting 
that the BAT data do generally support hard emission 
at higher energies, although at low signal-to-noise. Such 
hard emission could be due, on the other hand, to heav- 
ily obscured emission from a background AGN within 
the FOV. The spatial distribution of BAT emission is 
consistent with that from a point source in all bands. 

A 1651 (Fig. [6]): This cluster has a weak cool core, 
which means that while there is no significant temper- 
ature gradient in the c enter, the cooling tim e of the gas 
in the center is short (Hudson et al.l I2010D . Note that 
in a bimodal classification A1651 would not be consid- 
ered to have a cool core given its hig h central entropy 
of 90 keV cm 2 (jCavagnolo et al.ll2009t) . Given the simi- 
larity between its BAT data and that of A2029, an ob- 
scured AGN of similar flux could be responsible for the 
marginally detected positive flux in the higher energy 
bands. However, in this case the T+IC model is a sig- 
nificantly better fit than is the 2T model; A^ 2 improves 
by 9 (2-195 keV) and 5 (3-195 keV) over the IT and 
2T models. If there were no hard excess, the probabil- 
ity that the 6 highest energy bands measure flux above 
the thermal component, given that BAT fluctuations are 

Gaussian, is (|) , or 1.6%, which is not impressive in a 
sample of 59 clusters. The BAT spectrum is certainly 
suggestive, but considering the excess is not significant 
at the 3cr level for the 3-195 keV fit, and only just at 
this level in the 2-195 keV fit - without including sys- 
tematic uncertainties - we cannot claim to have detected 
a non-thermal component in this cluster. However, the 
evidence is perhaps strongest in this case, which is con- 
trary to the expectation that such an excess is most likely 
in a merging cluster, particularly one with a radio halo 
or relic. 

A2142 (Fig. [7]): As the hottest cluster in the sam- 
ple, the BAT is easily able to detect this cluster's 
high energy emission, which we might expect to ex- 
hibit a non-thermal excess sinc e it also hosts a radio 
halo (|Giovannini fc Ferettil [20001 . Both the T+IC and 
2T models indicate that hard excess emission may be 
present; in the latter case, the second temperature com- 
ponent is unphysically high, acquirin g the highest al- 
lowed temperature value. However, iNevalainen et all 
(2004) estimate that 2 Seyfert galaxy nuclei within 17' 
of the cluster center contribute ^30% of the hard band 
emission detected by Beppo-SAX; a similar amount of 
contamination would be expected in the BAT spectrum. 
Unfortunately, the XMM-Newton observation places this 
cluster right on the edge of the FOV, so over half (55.6%, 
based on a comparison with a pointed ROSATPSPC im- 
age) of the soft band emission is missing from the EPIC 
spectra. We rescale the XMM-Newton spectra to correct 
for the lost flux; the BAT source is equivalent to a point 
source, so it is not possible to correct the BAT emis- 
sion for the XMM-Newton FOV. The correction to the 
XMM-Newton flux could be off by a sizable factor if the 
E > 2 keV emission is distributed differently than the 



E < 2 keV emission where ROS AT is sensitive. The sig- 
nificance of the non-thermal excess here is only at the 2a 
level, mainly due to the poor statistics at XMM-Newton 
energies. While inconclusive, the BAT spectrum war- 
rants further analysis using better data below 12 keV. 

A3 112 (Fig El): Using both Cha ndra and XMM- 
Newton data, fBonamc nte et al.l (|2007|) have claimed to 
see both a hard and soft excess that is consistent with a 
non-thermal origin. If this is the correct interpretation of 
these spectra, the IC excess would be clearly detectable 
in the BAT spectrum given our sensitivity. While a non- 
thermal component is detecte d in our joint fits , it ha s 
well below the predicted flux of iBonamente et "all (2007); 
our 3<7 upper limit on the non-thermal normalization, us- 
ing a photon index T = 1.8 that matches their best- fit 
value, is 3 times lower than their estimate. The quality 
of our IT model fits is significantly less than for either 
the 2T or T+IC models; while those fits are of similar 
quality, the 2T fit yields physically reasonable temper- 
atures and lower x 2 values (Ax 2 ~ 3) than the T+IC 
model over both energy ranges. A non-thermal excess 
may in fact exist in this cluster, but a perhaps more likely 
scenario is that the ICM here is less isothermal than is 
typical in clusters, requiring several temperature compo- 
nents to adequately explain the cluster emission . The 
analysis of the Chandra data by iTakizawa et al.l ([2003D 
in fact demonstrates the multi-temperature structure of 
this cluster, which may be exaggerated by significant gas 
cooling outside the core. In any case, the BAT data 
do not argue strongly in favor of an IC interpretation 
for the excess emission above ~ 7 keV observed in the 
XMM-Newton data; as can be seen in FigureEl the power 
law component nearly ubiquitously overpredicts fluxes in 
the BAT spectrum. A more detailed exploration of the 
spatial and thermal structure at E < 12 keV is certainly 
warranted. 

A1367 ( Fig. EH): This cluster hosts a radio relic in its 
outskirts (|Gavazzi fc Trinchieri|[T98l . and so IC emis- 
sion is expected at some level in the radio relic re- 
gion; however, the XMM-Newton/ Swift extraction re- 
gion does not contain the relic, so we are unable to 
address the magnetic field s trength. Using RXTE, 
Hcnrikscn & Mushotzky] (|2001f) potentially detect a non- 
thermal component, although a two temperature fit bet- 
ter describes their spectrum. The marginally detected 
IC emission we see is consistent with their non-thermal 
flux, whether we use a photon index of 2.0 or their value 
(based on the spectrum of the radio relic) of 2.9. Our 
2T model fit, in the 2-195 keV band, is as good as the 
T+IC model fit, and given the marginally detected fluxes 
in the BAT bands, a 2T description of the ICM in this 
early stage, forming cluster cannot be ruled out. How- 
ever, the positive BAT fluxes and the consistency of our 
non-thermal fit with the analysis of the RXTE spectrum 
warrants future investigation of this cluster's hard X-ray 
emission. 

A2589 & Fornax: Neither of the BAT spectra of these 
clusters show particular evidence that that they have 
detected emission of any kind in any band. The first 
2 bands of A2589's spectrum are just inconsistent with 
zero flux at the 1<7 level, but a marginal detection in these 
bands is consistent with the thermal component. In both 
cases, the BAT spectrum is not sensitive enough to ex- 
clude the non-thermal component driven by the XMM- 
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Newton data; since the BAT data do not further con- 
strain the non-thermal component in these cases, we will 
not discuss these clusters further. 

4.3. Upper Limits 

While some evidence for non-thermal emission is 
present in several of the HIFLUGCS clusters, in none 
of these cases is a significant excess indicated by both 
the BAT and EPIC spectra that could not plausibly be 
explained by a multi-temperature state of the ICM. In 
many cases, the BAT spectra simply lacked the signal- 
to-noise to meaningfully constrain the existence of ex- 
cess emission; we therefore derive upper limits for a non- 
thermal component in our joint spectra. Three limits are 
presented for each energy range (2-195 keV and 3-195 
keV) considered: a 90% confidence level limit includ- 
ing systematic uncertainties in /qn and the EPIC back- 
grounds, as described in Section [2~T1 and two 3a limits, 
without systematic uncertainties included, for our fidu- 
cial photon index of T — 2 and for the best-fit value 
of r. After fitting for T, it is then fixed at that value 
when the upper limit is computed. The syste matic terms 
are in cluded in the 90% limits as described in lWik et al.l 
( 2009). Upper limits are reported as 20-80 keV fluxes in 
units of 10 -12 erg cm~ 2 s _1 in Tabled Note that when T 
is much steeper than 2, the power law component is con- 
strained only by the low energy spectrum and the 20-80 
keV flux limits are not reflective of the sensitivity of the 
BAT survey. In some instances, usually for lower tem- 
perature clusters, the 90% limit exceeds the 3cr limits; 
in these fits, the systematic uncertainties in /cn and/or 
the EPIC background dominates over the statistical un- 
certainty in the spectra. For example, in a low tempera- 
ture cluster lowering the EPIC backgrounds significantly 
hardens the spectra, while modifying /cn such that al- 
ready poorly constraining BAT fluxes are 10% higher, 
will allow a much larger IC-like component to fit the 
data than would be allowed statistically. In hotter clus- 
ters, adjusting the background has less of an effect on 
their spectral shape, and because they are hot they tend 
to be more significantly detected by the BAT, so that 
modifying /cn cannot drastically affect the non-thermal 
component. 

5. JOINT FITS TO STACKED EPIC-BAT SPECTRA 

In some clusters, as noted above, hints of a non-thermal 
excess are present, even if we cannot argue for their defi- 
nite detection. If the excess does exist in several clusters, 
but just below the detection threshold, we may be able to 
increase the signal-to-noise enough for a statistical detec- 
tion by stacking the cluster spectra. For simplicity, we 
stack only the EPIC-pn XMM-Newton spectra, which 
have the highest sensitivity especially at higher energies. 
Stacking the MOS spectra would be complicated by the 
variable pn/MOS cross-calibration factor and the fact 
that 3 of the cluster MOS spectra have been excluded 
from our analysis. Both the pn and BAT spectra are 
straightforwardly summed, as are the pn backgrounds, 
and their errors are propagated. Because the same re- 
sponse matrix is used for all the BAT spectra, we are 
able to use this unmodified file with the stacked spec- 
trum. To create an average response matrix for use with 
the stacked pn spectrum, we first multiply the individual 



redistribution matrices by their respective auxiliary re- 
sponse files, which contain the effective area per incoming 
photon energy. Then, a weighted average is performed 
on the new response files, with weighting factors pro- 
portional to each spectrum's 2-7 keV count rate. This 
procedure ensures that the final response matrix will best 
represent the instrumental response for the majority of 
photons. In any case, an unweighted response file was 
also created and no significantly different results were 
produced when using it. The CXB model normalizations 
were summed and included in the spectral fits. 

In all, we create 8 stacked spectra based on different 
groupings of the 59 HIFLUGCS clusters for which we 
have XMM-Newton data: "All" clusters, "Hot" (kT > 7 
keV, from the 2-12 keV fits), "Cool" (kT < 7 keV), 
"Radio" clusters hosting either a large-scale radio halo 
and/or relic or a smaller, central mini-halo, "No Radio" 
clusters that do not host any of these types of diffuse 
radio emission, non-cool-core clusters ("NCC"), strong 
cool core clusters ("SC C"), and weak c ool co re clusters 
("WCC"), as defined bv lHudson et~aLl (f2010h and listed 
in Table Q] These categories are designed to separate 
the sample into subgroups which might have different 
average levels of non-thermal emission. For example, IC 
emission must exist at some level in clusters with a radio 
halo or relic, but may not be present in clusters more 
generally. Thus, we might expect the "Radio" clusters 
to preferentially have non-thermal excesses, which are 
enhanced when they are stacked together and not diluted 
by the additional spectra from "No Radio" clusters that 
have no such excess. 

Because these clusters span a large range of temper- 
atures and redshifts, it is not appropriate to model the 
summed spectra with a single or even several temper- 
ature model for the thermal component. Instead, we 
build multi-temperature models from the previous spec- 
tral fits, for which we keep the spectral shape fixed and 
only allow the overall normalization to vary during fits 
to the stacked spectra. We consider the XMM-Newton- 
only single temperature fits (Table [5]) derived from 2- 
12 keV (1T X ,> 2 ) and from 3-12 keV (1T X ,> 3 ), and the 
single (lTj) and double (2Tj) temperature fits derived 
from the 2-195 keV joint spectra (Table [5]). To search for 
non-thermal emission in the stacked spectra, a power law 
model is added to represent the IC component and the 
normalization of the thermal model is allowed to vary. 
Ideally, the shape of the thermal component would be 
able to adjust to accommodate the IC signal, as it effec- 
tively does in the individual joint fits via the temperature 
parameter. However, the non-thermal flux below 12 keV 
will be small and should not cause the temperature to 
change in any significant way. For the 2Tj model, we 
want to avoid including unphysical temperature compo- 
nents that may have been driven by calibration features 
at the edges of the spectral range in the individual 2T 
fits. A low temperature (< 2 keV) component's emis- 
sion measure may cause < 2 keV emission to be signifi- 
cantly overestimated in order to better fit the gold edge, 
for example. Similarly, a slight under-subtraction of the 
XMM-Newton background or positive fluxes in the higher 
energy BAT bands may lead to unrealistically high tem- 
peratures. In Figure [TUl we plot the temperature values 
for this model relative to the lTj model temperatures. 
We have removed unphysical temperature components 
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TABLE 3 

Upper. Limits to 20-80 keV non-thermal Flux from EPIC and BAT Joint Fits 



Spectral Band: 2-195 kcV Spectral Band: 3-195 kcV 



Name 


90% a 


3cp = 2 a 




r 


3o- a 




90% a 


3 crp — 9 1 




r 


3a a 




A0085 


0.729 


6.66 


2 


.00 


6.66 




0.729 


6.66 


1 


.05 


5.20 




A0119 


1.689 


6.58 


2 


.12 


6.43 




1.689 


6.58 


9 


.14 


3.1x10" 


7 


A0133 


3.759 


6.88 


2 


.17 


6.40 




3.759 


6.88 


2 


.14 


4.71 




NGC507 


1.725 


1.77 


1 


.90 


2.13 




1.725 


1.77 


1 


.97 


2.08 




A0262 


0.804 


4.05 


2 


.00 


4.05 




0.804 


4.05 


3 


.77 


0.15 




A0400 


0.645 


2.26 


9 


.02 


2.1x10" 


8 


0.645 


2.26 


8. 


.80 


6.9x10- 


7 


A0399 


1.818 


4.04 


2 


.07 


3.62 




0.032 


4.04 


4 


.53 


0.01 




A3112 


5.420 


8.98 


2 


.00 


8.98 




5.420 


8.98 


2 


.18 


6.98 




Fornax 


3.365 


3.13 


2 


.00 


3.13 




3.365 


3.13 


2. 


.18 


2.67 




2A0335 


1.626 


6.80 


2 


.00 


6.80 




1.626 


6.80 


2 


.18 


4.31 




IIIZw54 


2.626 


4.32 


2 


.14 


3.48 




2.626 


4.32 


2 


.33 


2.63 




A3 158 


2.453 


5.08 


1 


.93 


5.24 




2.453 


5.08 


1 


.98 


4.56 




NGC1550 


2.478 


1.66 


2 


.06 


1.94 




2.478 


1.66 


2 


.39 


1.60 




EXO0422 


4.171 


5.18 


2 


.01 


5.14 




4.171 


5.18 


2 


.03 


4.98 




A3266 


5.131 


7.25 


1 


.98 


7.43 




5.131 


7.25 


1 


.89 


7.87 




A0496 


0.733 


6.28 


2 


.00 


6.28 




0.733 


6.28 


2. 


.11 


3.33 




A3376 


5.582 


6.01 


2 


10 


5.03 




5.582 


6.01 


2 


.02 


5.77 




A3391 


4.609 


6.30 


2 


.00 


6.25 




4.609 


6.30 


2 


.00 


6.08 




A3395s 


3.548 


3.52 


2 


.00 


3.44 




3.548 


3.52 


2. 


.37 


2.48 




R1504 


9.934 


5.87 


2 


.01 


5.72 




5.826 


5.87 


2 


.62 


31.90 




A0576 


3.671 


6.44 


2 


.00 


6.44 




3.671 


6.44 


2 


.02 


6.36 




A0754 


4.332 


7.61 


2 


.00 


7.61 




4.332 


7.61 


1 


.98 


13.77 




HydraA 


4.169 


8.72 


2 


.08 


8.55 




4.169 


8.72 


2 


.06 


5.78 




A1060 


0.469 


7.30 


2 


.19 


7.11 




0.469 


7.30 


2 


.61 


1.49 




A1367 


5.676 


8.24 


2 


.00 


8.37 




5.676 


8.24 


2 


.00 


7.30 




MKW4 


3.453 


1.98 


1 


.93 


2.07 




3.453 


1.98 


1 


.69 


3.80 




ZwC11215 


1.585 


5.45 


2 


.00 


7.27 




1.585 


5.45 


2 


.26 


2.84 




NGC4636 


0.886 


0.61 


2 


.00 


0.61 




0.886 


0.61 


2. 


.16 


0.63 




A3526 


3.709 


11.23 


2 


.00 


11.59 




3.709 


11.23 


2 


.16 


2.03 




A1644 


4.915 


8.65 


2 


.00 


8.65 




4.915 


8.65 


2 


.08 


7.53 




A1650 


4.306 


6.63 


2 


.00 


6.62 




4.306 


6.63 


2 


.00 


6.71 




A1651 


8.497 


10.82 


2 


.00 


10.83 




8.497 


10.82 


1 


.96 


11.06 




Coma 


1.509 


6.55 


2 


.03 


6.48 




1.509 


6.55 





.47 


3.20 




NGC5044 


1.831 


1.92 


2 


.14 


1.46 




1.831 


1.92 


2 


.45 


0.89 




A3558 


0.229 


9.20 


2 


.53 


1.85 




0.223 


9.20 


2 


.58 


0.64 




A3562 


2.387 


5.52 


4 


.98 


1.3x10" 


3 


2.311 


5.52 


8 


.97 


7.8x10- 


7 


A3571 


0.442 


8.52 


2 


.39 


3.51 




0.442 


8.52 


9 


.46 


1.3x10" 


7 


A1795 


2.143 


7.22 


2 


.68 


1.09 




2.143 


7.22 


2 


.10 


2.99 




A3581 


2.351 


2.76 


2 


.00 


2.75 




2.351 


2.76 


2 


.34 


1.70 




MKW8 


3.508 


4.99 


2 


.04 


4.80 




3.508 


4.99 


2 


.05 


4.92 




A2029 


7.633 


13.33 


2 


.00 


13.31 




7.633 


13.33 


1 


.82 


11.50 




A2052 


3.430 


6.95 


2 


.00 


6.95 




3.430 


6.95 


2 


.14 


4.60 




MKW3S 


3.280 


5.96 


2 


.00 


5.96 




3.280 


5.96 


1 


.95 


5.34 




A2065 


3.842 


5.59 


2 


.00 


5.60 




3.842 


5.59 


2 


.00 


5.95 




A2063 


2.573 


6.58 


2 


.13 


5.81 




2.580 


6.58 


2 


.16 


5.01 




A2142 


11.494 


13.51 


1 


.95 


13.48 




11.494 


13.51 


1 


.98 


14.65 




A2147 


2.203 


5.50 


2 


.01 


5.52 




2.203 


5.50 


8. 


.62 


2.1x10" 


6 


A2199 


3.768 


7.53 


2 


.20 


8.05 




3.768 


7.53 


2 


.17 


5.89 




A2204 


5.809 


11.75 


2 


.09 


9.30 




5.809 


11.75 


2 


.05 


8.11 




A2256 


3.013 


7.38 


2 


.28 


5.07 




3.013 


7.38 


2 


.36 


2.98 




A2255 


1.321 


2.73 


6 


.24 


2.3x10" 


5 


1.321 


2.73 


9 


.41 


2.5x10" 


7 


A3667 


1.932 


12.00 


2 


.39 


4.01 




1.932 


12.00 


2 


.20 


2.79 




S1101 


1.009 


2.23 


2 


.43 


1.65 




1.009 


2.23 


2 


.08 


2.13 




A2589 


5.045 


5.82 


2 


.00 


5.82 




5.045 


5.82 


2. 


.02 


6.42 




A2597 


2.620 


6.10 


2 


.00 


6.09 




2.620 


6.10 


2 


.41 


2.90 




A2634 


1.806 


4.34 


7 


.28 


3.3x10- 


6 


1.806 


4.34 


9 


.50 


4.7x10- 


7 


A2657 


0.613 


4.58 


2 


.00 


4.58 




0.613 


4.58 


7. 


.28 


2.6x10- 


5 


A4038 


3.117 


7.43 


2 


.00 


7.43 




3.117 


7.43 


2 


.12 


4.73 




A4059 


1.981 


3.82 


2 


.44 


2.44 




0.060 


3.82 


9 


.85 


6.1x10" 


8 


a 20-80 keV, 


lO" 12 


erg cm~ 2 
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from both the 2Tj model; the best-fit single tempera- 
ture model is used in place of the 2T model for those 
clusters, which are represented by blue circles in Fig- 
ure [101 Unphysical temperature components were found 
to have kT > 16 keV and kT < 2.1 kcV, if their lTj tem- 
perature is greater than 3.5 keV. In general, this latter 
cut eliminates temperature components that significantly 
over-predict the 0.5 keV < E < 2 keV emission. 

Thermal and thermal plus non-thermal fits to the 
stacked spectra are given in Table |U Considering only 
the fits to data with E > 3 keV, which excludes the most 
problematic region of the spectra, we find no evidence at 
the statistical 90% level for a non-thermal component in 
any of the stacked spectra except in the case of "Radio" 
clusters. In the table, the normalization of the thermal 
component in the "TModei-only" fits is not shown, only its 
\ 2 value for comparison purposes. For the "TModei+IC" 
fits, the photon index is fixed to T = 2 as was done pre- 
viously for the joint fits. The last 3 columns report the 
"TModei+IC" fits with T as a free parameter; however, 
its value is fixed when errors are computed. In this case, 
the photon index was initialized as T = 2, so for spec- 
tra with no particularly strong indication of non-thermal 
emission, the best-fit normalization was set to zero and 
the photon index kept at or near its initialized value; this 
explains why so many of the "best-fit" photon indices 
presented in the table are '2.00.' In the case of large val- 
ues of r > 3, the non-thermal component is attempting 
to either represent incompletely modeled soft emission 
from low temperature gas or correct an imperfectly cali- 
brated gold edge. Even though these normalizations are 
large and quite significant, they are so steep that the 
flux at hard energies is negligible and does not represent 
an IC excess. If < 2 keV emission were included in the 
fits, these large T values would disappear as they would 
vastly over-predict the soft emission. 

In Figure [TTJ the joint fit for the stacked spectra of all 
59 clusters is shown with the lTx,>2 model. The best- 
fit model normalization agrees with its expected value to 
better than 1%, as do all the model fits without an IC 
component, indicating that the average pn response is ac- 
curate. Also, a difference in spectral shape appears below 
3 keV, visible in the residuals, that highlights the prob- 
lem with including this emission in the fits. The BAT 
data are well represented by this model, even though the 
temperature models were derived from fits to the XMM- 
Newton spectra alone. The regular pattern in the BAT 
residuals is likely real, and is apparent in most of the 
spect ra of hot clusters such as Coma (see W ik et al.l 
1201 If ). When considering only one cluster, it seemed rea- 
sonable that this residual pattern could simply be due to 
chance. The pattern reappears in many of the individ- 
ual joint fits however, indicative of a systematic problem. 
Because the BAT flux calibration is dominated by nor- 
malizing to the Crab flux in each band, these fluxes are 
really only accurate for objects with a spectral slope sim- 
ilar to the Crab's. At these energies, cluster spectra are 
quite steep even for the hottest temperatures, so some 
miscalibration would be expected. Most likely, the first 
and possibly second energy bands have underestimated 
fluxes, owing to the rapid rise of the instrumental re- 
sponse with energy; clusters have proportionately more 
emission at the lower energy part of the band than does 
the Crab, and so the internal band response is miscali- 
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kT (keV) [1T Model (2-195 keV)] 

Fig. 10.— The two temperature values in the 2T fits (2-195 keV) 
with respect to the single temperature fit for each cluster over the 
same range. High (triangles, black) and low (circles, red) temper- 
ature values in the 2T model for the same cluster are connected by 
dotted lines for clarity. Blue circles are IT fits to clusters with an 
unphysical best-fit 2T model, with either unrealistically high tem- 
peratures (kT > 16 keV) or a low temperature component that 
over-predicts the emission below 2 keV. (These have kT < 2.1 keV 
for clusters with single temperatures of 3.5 keV or hotter.) These 
excluded regions are indicated by dashed lines in the figure, and 
the diagonal dashed line represents equality of IT and 2T temper- 
atures. These temperature values are used to build the 2Tj model 
used in fits to the stacked spectra. 




Energy (keV) 

Fig. 11. — The stacked spectrum of all 59 clusters with the com- 
bined single temperature model fit (lTx,>2)- The EPIC-pn spec- 
trum (E < 12 keV) and BAT spectrum (E > 14 keV) are shown in 
the top panel, and their residuals in the lower panel. The CXB con- 
tribution appears below the EPIC-pn data. The problems between 
2-3 kcV (described in the text) clearly show up in the residuals, as 
does a potential problem with low energy BAT fluxes. The com- 
bined single temperature model determined from the 2—12 keV fits 
is sufficient to explain the summed BAT spectrum; no non-thermal 
excess is obvious. 

brated - weighting the higher energy part of the response 
more strongly than is appropriate for thermal emission. 
While this certainly affects our results, the only solu- 
tion is to develop a detailed response matrix model for 
the survey data. Unfortunately, the detailed spectral 
response for the Swift survey data currently has much 
larger uncertainties than the Crab spectrum itself. 

In general, the addition of a non-thermal component 
to these spectra does not significantly improve the fits 
in Table 0] except for the "Radio" - and to a lesser de- 
gree the "Hot" - subsamples. The "All," "Cool," "No 
Radio," "WCC," and "SCC" stacks are found to lack 
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a physically plausible (r < 3) non-thermal component 
at the statistical-only 90% level. For the "Hot" and 
"NCC" sample fits, the IC component improves fits us- 
ing the lTx,>2 and lTj thermal models, but not the 
lTx,>3 and 2Tj models. In contrast, the "Radio" clus- 
ter sample T+IC fits are not only clearly improved over 
the thermal-only fits, but the IC component is significant 
at the 90% level regardless of the thermal model consid- 
ered. Also, except in the case of the 2Tj model, allowing 
the IC index to vary shows that the slope is consistent 
with expected indices (2-2.5). The exact value of the in- 
dex favored by the data should not be taken to represent 
the true shape of the non-thermal component, however, 
since it is most strongly influenced by the 2 keV < E < 3 
keV data and thus biased by gold edge calibration issues 
and incompletely modeled cool gas. 

The best-fit non-thermal plus Tx,>2 model for the Ra- 
dio clusters is shown in Figure 1121 For comparison, the 
Tx,>2 and 2Tj fits with no IC component are shown in 
Figure [TBI and Figure Qj] respectively. The non-thermal 
component, plotted as a dotted line in the figure, be- 
comes competitive with the thermal emission in the 35- 
50 keV band, where a somewhat significant excess is 
present in thermal-only model fits. By contrast, the "No 
Radio" subsample shows no evidence for an excess at 
hard energies (Fig. [TS]). 

Ignoring systematic uncertainties, the non-thermal sig- 
nature is detected for the Radio clusters with 5.1(7 con- 
fidence using the Tx,>2 model and 2.4er with the Tx,>3 
model. Including a conservative /cn uncertainty of 10%, 
which assumes the average cross-calibration is incorrect 
by that much, reduces the significances to 3.3<r and 0.7tr, 
respectively. While the detection is quite robust with the 
Tx,>2 thermal model, both the normalization and signif- 
icance of the detection degrades when using the Tx,>3 
model. Since the temperatures making up the models 
are not very different (see Fig. 1101) . the primary driver 
for this change must come from removing the 2-3 keV 
data when fitting with the Tx,>3 model. The decreased 
significance should not necessarily be a concern as the 
lowest energies statistically dominate x 2 - However, the 
factor of 2 drop in flux associated with the IC compo- 
nent between the Tx,>2+IC and Tx,>3+IC fits suggests 
that the power law is not driven by broadband IC emis- 
sion but instead by features around 2-3 keV. We cannot 
therefore reasonably claim an ensemble detection of ex- 
cess hard emission at the 3<r level given the systematic 
and modeling uncertainties. 

6. IMPLICATIONS AND DISCUSSION 

In this work, we characterized the hard X-ray emis- 
sion from HIFLUGCS, a sample of the brightest galaxy 
clusters outside the Galactic plane. For the 59 out of 
64 clusters with usable XMM-Newton data, we searched 
for excesses over the thermal emission from gas in the 
ICM in data from the 58-month Swift BAT all-sky sur- 
vey. EPIC and BAT spectra were extracted from identi- 
cal regions and carefully calibrated to allow straightfor- 
ward joint fits that simultaneously constrain the thermal 
and non-thermal emission in both spectra. We first con- 
sidered fitting over an energy range of 2-195 keV but 
found that low temperature gas and the gold edge in 
the XMM-Newton spectra could lead to false detections. 
Ignoring the 2-3 keV data resolved this issue, although 




Energy (keV) 

Fig. 12. — The stacked spectrum of all clusters with large-scale, 
diffuse radio ha los or relics. General features of the plot are the 
same as Figure 1111 The dotted line represents the best-fit non- 
thermal model with photon index T = 2, and the upper solid line 
shows the lTx,>2 thermal model. 
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Fig. 13. — The stacked spectrum of all clusters with large-scale, 
diffuse radio halos or relics with the combined single temperature 
model fit (lT x >2)- General features of the plot are the same 
as Figure 1111 A slight excess is apparent in the BAT spe ctru m, 
due to either a non-thermal spectral component (see Fig. 112b or 
a signific ant multi-temperature structure in the individual clusters 
(see Fig. [14J 
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Fig. 14. — The stacked spectrum of all clusters with large-scale, 
diffuse radio halos or relics with the combined double tempera- 
ture mod el fi t (2Tj). General features of the plot are the same 
as Figure 1111 The combined 2Tj model can explain most of the 
slight excess seen when the single temperature model (lTx >2) is 
considered. 
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TABLE 4 

Fits to Stacked EPIC and BAT Spectra 



Sample 




TModol^only 


T Mode 


+IC 




Tlvlodcl+IC, rf ree 


(number) 


Tj^odcl 


X 2 /dof 


Norm. a 


X z /dof 


r 


Norm. a 


X^/dof 


All 


lTx,>2 


1218.34/1606 


< 0.0165 


1217.15/1605 


6.94 


1 0606+!!- 1 ? 5 ! 

uuu -0.1657 


1107.44/1605 


(59) 


1T X ,>3 


954.95/1406 


< 0.0041 


954.95/1405 


2.00 


< 0.0041 


954.95/1405 




ITj 


1218.62/1606 


< 0.0162 


1217.56/1605 


7.16 


1 9145+ - 1926 


1110.99/1605 




2Tj 


1225.78/1606 


< 0.0062 


1225.78/1605 


2.00 


< 0.0062 


1225.78/1605 


Hot 


lTx,>2 


1037.57/1606 


_ „„„„.(.n 004Q 

0.0093±°;°°g 


1027.64/1605 


2.21 


m47+ - 0067 
u.ui4<_ 0067 


1024.50/1605 


(12) 


lTx,>3 


862.44/1406 


< 0.0082 


861.68/1405 


0.83 


< 9.3xl0~ 5 


860.09/1405 




ITj 


inoT 1 -\ eric 

1037.76/1000 


n nncQ+0.0049 
U.UU8»_ 004g 


1028.96/ 1605 


2.22 


0.0140+°^7 


1025.95/1605 




2Tj 


1033.29/1606 


< 0.0093 


1031.20/1605 


8.02 


0.5347+° :2 °|| 


1015.51/1605 


Cool 


lTx,>2 


1101.49/1606 


< 0.0069 


1101.49/1605 


2.00 


< 0.0069 


1101.49/1605 


(47) 


lTx,>3 


892.09/1406 


< 0.0036 


892.09/1405 


2.00 


< 0.0036 


892.09/1405 




ITj 


1101.51/1606 


< 0.0078 


1101.51/1605 


7.62 


0.97481°-^^ 


1047.00/1605 




2T.j 


1098.84/1606 


< 0.0034 


1098.84/1605 


2.00 


< 0.0034 


1098.84/1605 


Radio 


lTx,>2 


1007.22/1605 


-i_n nn^^ 

0.0137±g;8gil 


990.54/1604 


2.29 


0265+ - 0080 

U.UZOD_q q 080 


977.31/1604 


(15) 


1T X ,>3 


835.82/1405 


U.UUD1_q 005g 


832.91/1404 


2.01 
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Fig. 15. — The stacked spectrum of all clusters without diffuse 
radio emission shown with the combined single temperature m odel 
fit (lTx,>2)- General features of the plot are the same as Figure lTTl 



a somewhat weaker constraint on the thermal compo- 
nent reduced our overall sensitivity. From the 3-195 keV 
fits, six clusters were found to have marginal evidence 
for a non-thermal excess, although none of these were 
deemed significant enough to claim a detection, espe- 
cially considering systematic uncertainties in the EPIC 
background and EPIC-BAT cross calibration normaliza- 



tions. We then stacked the spectra to look for a signifi- 
cant statistical detection of non-thermal emission in the 
HIFLUGCS sample. Unfortunately, the stacked spectra 
revealed no definitive excess. Stacking subsamples of the 
HIFLUGCS clusters returned similar results, except for a 
tantalizing but very marginal detection of a non-thermal 
component in the stacked spectrum of all clusters that 
host radio halos and/or relics or mini-halos — the very 
clusters that are most expected to have detectable IC 
emission. 

6.1. Comparison to Previous Studies 

The lack of definitive hard X-ray excesses in our 
individual clusters is consistent with the most recent 
searches with Suzaku, INTEGRAL, and Swift, though 
somewhat less so with those of RXTE and Beppo-SAX. 
Ignoring the Coma cluster, whose controversial hard 
energy emission is discussed at length elsewhere (e.g., 
iWik et all 1201 1[ ). our analysis is not clearly inconsis- 
tent with any previous observations, particularly given 
that the possible existence of low-level, extended non- 
therm al emission has n ot been considered in detail here 
(as in IWik et al1l20ll . which RXTE and Beppo-SAX 
in particular would be sensitive to given their large 
FOVs. For the c lusters in our sample al so observed 
by RXTE, A3667 (IRephaeli fc Gruberl 120041) and A2256 
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(|Rephaeli fc Gruberil200l . our upper limits agree with 
analyses of their data, at least considering the two- 
temperature interpretation allowed for A2256, regardless 
of the distribution of emission. The recent RXTE detec- 
tio n of non-th e rmal emission in NGC 5044 below 15 keV 
by iHenriksenl (|2011| ) lies below our detection threshold 
at higher energies. For several of the clusters observed 
with Beppo-SAX an d found to host non-therm al emis- 
sion, such as A2256 (iF usco-Femiano et al.ll2005D, A219 9 
(|Kaastraet al.lll999f ). and A3526 (|Molendi et atll2002h . 
our upper l imits fall below their measured inverse Comp- 
ton fluxes. iKaastra et al.l ilT9~99i claim an extended non- 
thermal halo for A2199 between 0.5 and 1.5 Mpc, which 
is not inconsistent with its larger size at high energies 
(14-20 keV, see Fig. [I]); however, due to the low S/N 
of the detection, this extent is also indistinguishable 
from that of a point source. Upper l imits from Suzaku 
for clusters A3667 jNakazawa et all I2009Q and A3376 
(|Kawano et al. 2009]) are obviously consistent with these 
results. 

Similar studies of c lusters detected by the BAT 
(|Ajello et al.ll2009Ll2010h have also failed to find definitive 
non-thermal ex cesses. The only dis crepa ncy is for A3667, 
for wh ich both lAjello et al.l (|2010t ) and iNakazaw a et al.l 
( 2009) detect high temperature (kT ~ 15 keV) gas near 
the center. While we do not see strong evidence for a sig- 
nificant high temperature component like this - although 
our 2T, 2-195 keV fit does suggest a significant amount 
of hot gas (kT ~ 9 keV) - the elongated shape caused by 
its ongoing merger requires a more detailed analysis to 
more accurately extract its BAT fluxes to properly assess 
this high temperature component. In any case, a note- 
worthy difference between the methodology here and in 
iWik et all (poTl with that of lAjello et all (120091 [2010h 
is our use of the technique developed bv Ikenaud et all 
(2006) to recover extended source fluxes from coded mask 
observations. This procedure allows for a more direct 
spatial comparison between soft and hard X-ray spectra 
such that no assumptions about the extent of hard band 
data need to be made; however, the low relative extent 
and signal-to-noise generally achieved makes this advan- 
tage critical only for the largest, brightest clusters such 
as Perseus and Coma. 

While some excesses in the stacked spectra are tan- 
talizing, equally good, and sometimes better, fits result 
when the 2Tj model is used. Since only the normaliza- 
tion is allowed to vary in these fits, it is hard to justify 
why the addition of an IC component really provides a 
better description of the data, especially if the improve- 
ment in x 2 is minor. Note that this comparison is only 
fair because the 2T models are all physically reasonable 
descriptions of the ICM, otherwise we may be inappropri- 
ately modeling non-thermal emission with an incorrect 
thermal component. The upper limits on non-thermal 
emission in the stacked spectra, when applied on aver- 
age to the clusters making up the stacked sample, are 
more constraining than limits from individual fits. The 
typical 90% confidence level upper limit on the cumula- 
tive IC flux in the stacked spectra is 10 -11 erg cm~ 2 s _1 
in the 20-80 keV band, which translates to an average 
maximum flux per cluster 4 to 8 times lower than this 
limit. 

These results are in conflict with an analysis of 
a similar sample of clusters observed by Beppo-SAX 



(Ncvalain en et all 120041 ) . which found systematic if 
marginal excesses for merging clusters. Actually, these 
previous IC flux estimates are not unlike our results in 
the 2-195 keV range, as are the temperatures of the ther- 
mal component for clusters in both our and their sam- 
ples. However, over the 3-195 keV energy range, the 
90% error interval for nearly all the excesses include zero. 
This result is at least partl y du e to slightly higher best-fit 
temperatures (see Section |3~T1 for a more detailed discus- 
sion). Since clusters are not isothermal, harder spectra 
such as those from the BAT will contain proportionately 
more photons from higher temperature gas. An exam- 
ple of this bias c a n be seen in the stacked spectrum of 
iNevalainen et all (|2004f) : they observe a highly signifi- 
cant non-thermal excess, but the steep IC component 
necessary to explain it would lead to detectable amounts 
of non-thermal emission at softer energies, which is not 
seen. The authors interpret this as evidence that the 
non-thermal emission is significantly extended. Our BAT 
data test this possibility, as extended emission is both de- 
tectable and not detected by the BAT beyond that pro- 
duced by the thermal gas. Thus, it is unlikely that the 
non-thermal emission is very highly extended and strong. 
In fact, the steep excess in their spectrum is exactly what 
would be expected for a strongly multi-temperature ther- 
mal structure - which naturally results when many clus- 
ters spanning a broad range in temperature are summed 
- that is modeled as a single temperature component, 
which is what they do. When we model our stacked BAT 
spectrum this way, we find a temperature consistent with 
the average temperature of our clusters and a very sig- 
nificant, steep (r ~ 2.8) power law component, identical 
to their best-fit photon index. But the thermal compo- 
nent, determined at hard energies, will be more highly 
weighted by hotter clusters, whose emission dominates. 
If a single temperature component is used to model the 
thermal emission for such a summed spectrum, then at 
the very least the temperature needs to be fixed to the 
weighted-average value in the band in which the hard 
excess emission is expected to be found. For example, 
in our sample, the count rate weighted-average temper- 
ature jumps from 5.6 keV, when weighted by the 2—7 
keV count rate, to 7.1 keV when the 14-50 keV count 
rate is used. Here we have employed the temperatures 
determined from the 2-12 keV fits. Even so, the ex- 
act value of the temperature is less important than the 
fact that the highly multi-temperature composite spec- 
trum no longer looks like a single temperature plasma. 
The proper procedure is to use a truly multi-temperature 
model based on the temperatures of the constituent clus- 
ters, as we have done. We suspect that, if the thermal 
co mponent is similarly mo deled for the stacked spectrum 
of Ncvalai nen et al.l ( 2004f ). the non-thermal excess will 
be reduced; however, it is unlikely that all of their excess 
would disappear. 

6.2. Implied Magnetic Field Strengths 

Our most suggestive result from the various stacked 
subsamples, that clusters hosting a radio halo or relic 
have the the most significant indication of a non-thermal 
excess on average, is also the least surprising of possi- 
ble outcomes. Because radio halos and relics are as- 
sociated with mergers, which also produce shocks and 
multi-temperature gas distributions, the more appropri- 
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ate thermal model to use might be the 2Tj model, al- 
though even in this case a non-thermal component im- 
proves the fit. Assuming the power law component rep- 
resents true IC emission, we can determine the average 
B field in these clusters from the IC and synchrotron 
flux ratio as outlined in iWik et all (|2009f ). Unfortu- 
nately, not all of the diffuse radio emission of these clus- 
ters falls within the FOVs of the observations, so the 
measured IC normalization in Table [4] is not the best 
value we can use. We therefore sum a new subsam- 
ple of cluster emission, excluding those with radio emis- 
sion outside their XMM extraction regions, most no- 
tably the relics in A3667, A1367, and A3376. For the 
radio halo/relic clusters, we are l eft with radio emis- 
sion f rom Coma (4.4 Jy at 74 MHz (lHanisch fc Ericksonl 
1980), vign etted for the extrac tion region), A754 (4 Jy 
at 74 MHz dKassim et alJl20qil) ). A3562 (220 mJy at 240 
MHz (iGiacintucci et al.ll2005D ). A2256 (100 Jy at 22 MHz 
(|Costain et al.lll972h . halo and relic emission combined) , 
A2255 (475 mJy at 150 MHz (jPizzo & de Bruvnl [20091) . 
halo emission only) , and A0399 (16 mJy at 1400 MHz 
(jMurgia et al.l 120101 ) ) . Extrapolating these flux densi- 
ties to a common 74 MHz given their individual spectral 
indices, we find a total flux of 32 Jy. Low frequency 
flux densities have been used where available to mitigate 
the effect of spectral curvature on the IC/synchrotron 
flux ratios. For the mini-halos, low frequency measure- 
ments are less common: A 2142 (18 mJy at 1400 MHz 
(iGiovannini fc Ferett.il 120001) V R1504 (121 mJy at 327 
MHz dGiacintucci et al.ll2011h h A2029 (18.8 mJy at 1400 
MHz (iMurgia et all 120091) ). 2A0335 (22.6 mJy at 1500 
MHz (Sar azin et al.l 1995D). The two remaining mini- 
halos, A2204 ([Sanders et all I2009D and A2052 (Clarke, 
priv. comm.), do not have published diffuse fluxes. As- 
suming spectral indices of a = 1.3 {F u ex v~ a ) and an 
average flux density of 18 mJy at 1400 MHz for all 6 
clusters, a total flux of ~ 5 Jy at 74 MHz is found. Com- 
bining all these flux densities and using the power law 
normalization range in Table [4] (20-80 keV fluxes of (1.3- 
3) x 10~ n erg s _1 cm~ 2 ) yields an average B ~ 0.13- 
0.19/xG. 

Because mini-halos are confined to the cool cores of 
more relaxed clusters, completely unlike large-scale ha- 
los and relics, it may be more appropriate to consider 
the two classes separately. Surprisingly, the fit to the 
new radio halo/relic stacked spectrum, made up of the 
clusters listed above, completely disfavors the addition of 
a non-thermal component. Taking the 90% statistical- 
only upper limit to the IC flux of 5.7 x 10~ 12 erg s _1 
cm -2 , we estimate B > 0.26/xG on average in the cen- 
tral regions of these clusters. In contrast, the best-fit IC 
flux to the stacked mini-halo spectrum is 1.9 x 10 -11 erg 
s _1 cm -2 , significant at nearly the 4cr level (statistical- 
only), implying a magnetic field strength of 0.08/xG on 
average. For both cases the lTj thermal model is used 
over the 2-195 keV band. Individually, the best-fit IC 
fluxes of A3562, A0399, and A2255 are quite low (im- 
plying B > 10/iG) with 90% (stat.+sys.) lower limits 
of 0.06 fiG, 0.04 /xG, and 0.09 /xG, respectively. The 
remaining halo/relic clusters are more accommodating 
of a non-thermal component, with best-fit /lower limit 
B field strengths of 0.21/0.11 fiG (A0754), 0.54/0.21 
fiG (Coma), and 0.46/0.33 x*G (A2256). To compute 



these and the following B field strengths, we use the 
average of the best-fit normalizations in Tables [5] and 
[6] and the larger upper limit from Table (3[ The mini- 
halo spectra typically favor the additional power law 
component, leading to nominal best-fit IC fluxes not 
much lower than their upper limits and corresponding B 
field strengths of 0.16/0.12 fiG (2A0355), 0.06/0.04 /xG 
(A2142), 0.08/0.06 /xG (A2204), 0.06/0.05 xxG (A2029), 
0.08/0.07 /xG (A2052), and 0.10/0.05 /xG (R1504), as- 
suming a spectral index of 1.3 and radio flux densities of 
- 18 mJy at 1.4 GHz for A2052 and A2204. 

These results suggest that mini-halos may involve lower 
B fields, more easily allowing the detection of IC emission 
as compared with larger halos and relics. This general 
conclusion is consistent with the measurement of non- 
thermal emission with INTEGRAL associated with the 
mini-halo in the core of the Ophiuchus cluste r that im- 
plies B = 0.05-0.1/xG (jNevalainen et alj |2009). While it 
is tempting to associate the IC flux with the entire mini- 
halo sample, the vast majority of the BAT flux originates 
with A2029 and A2142 (three-quarters of the 14-24 keV 
emission and > 90% of it at the higher energies). These 
two cases are discussed individually in Section l4~2l A2142 
is particularly problematic, in that the XMM-Newton ex- 
traction region only contains half the cluster, which adds 
significant systematic uncertainty to the cross-calibration 
with the BAT spectrum, and in that the classification of 
the diffuse radio emission as a mini-halo, as opposed to 
radio galaxy jets, is not ironclad. We are therefore pre- 
cluded from drawing strong conclusions from this result, 
which is robustly driven by only a single cluster: A2029. 
In any case, the non-thermal component in the mini-halo 
spectrum is consistent with and somewhat driven by the 
BAT measurements, suggesting that mini-halos may be 
better targets of future hard X-ray observations, espe- 
cially given their smaller size, than traditional halos and 
relics. 

It may not be surprising that IC emission was not de- 
tected definitively in the radio halo/relic clusters; direct 
measurements of cluster magnetic fields through Faraday 
rotation measure (RM) studies typicall y find line-of-sight 
B fie lds on the order of several /xG (jGovoni k. Ferettl 
120041) . Similar high values of B are sug gested by the 
stabi lity of cold fronts in merging clusters ([Keshet et al.l 
I2010I ). although the flow may locally amplify the fields 
in these regions, so they would not be representative 
of the global average field strength, even in mini-halos. 
Also, RM magnetic field strengths could be biased high if 
stronger fields are correlated with denser gas, since RM 
observations are really measuring the electro n density- 
weigh ted value of B along the line of sight dPetrosianl 
120011 ). Such explanations, while entirely reasonable, were 
primarily developed to explain the lower values of B im- 
plied by earlier IC detections, some of which have been 
more recently c alled into question (e.g., w ith Suzaku, 
iNakazawa etaL|[2007l [20091 : 1 Wik et al J 120091) . However, 
if our low significance evidence for IC emission associated 
with mini-halos can be corroborated, then such low val- 
ues of B may in fact be common. Our current sensitivity 
to IC emission with either pointed or survey observations 
can only detect non-thermal emission in clusters with ra- 
dio halos if the magnetic fields are <0.2 /xG. Note that it 
is possible to observe much fainter IC emission at lower 
X-ray energies, and thus measure larger B fields, in ra- 
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dio relics that are significantly displaced from the X-ra y 
luminous gas in cluster centers ()Finoguenov et aLll2010f) . 

6.3. Dynamical Importance of the Non-thermal 
Electron Population 

Given that we have searched for and provided lim- 
its on IC emission from relativistic electron populations 
in a complete sample of clusters, what does that tell 
us about the general energetic importance of this pop- 
ulation relative to the thermal electron population of 
clusters in the low redshift universe? Hydrostatic clus- 
ter mass estimates - employed by studies using clusters 
as cosmological probes - typically assume a negligible 
amount of non-thermal pressure in central regions. In 
the following, we attempt to confirm this presumption. 
A power law IC spectrum is thought to originate from 
an energy distribution of relativistic electrons propor- 
tional to 7" p , where 7 is the electron Lorentz factor and 
p = 2a + 1 = 2T — 1 . The total energy in these electrons 
depends on the range of 7s and is particularly sensitive 
to the lower cut off energy where mo st of the electrons re- 
side. Following lMurgia et all (|2010f ). a rare case in which 
rough observational estimates exist, we take 7 m i n ~ 300 
and 7 m ax ~ 3 x 10 4 . The large uncertainty in these 
choices, which are impossible for us to estimate since 2— 
200 keV energies correspond to 1500 < 7 < 1.5 x 10 4 , 
limits the following discussion to order of magnitude es- 
timates at best. 

Our goal is to compare the typical non-thermal pres- 
sure, Pnt = unt/3, where «nt is the energy density 
of the non-thermal population, to the typical thermal 
pressure inside our extraction regions. For simplicity, 
we ad opt the universal pressure profile of lArnaud et al.l 
(2010) from which the thermal pressure inside the region 
(taken to be ~ | of the extraction radius) can be derived 
from the mass estimates at i?5oo of Zha ng et al.l (|2011h . 
In individual clusters, the 90% upper limits on the IC 
flux typically limit the non-thermal electron to thermal 
electron pressure to less than a few tens of percent. The 
limits are stricter for massive (hot) clusters (< 10%), pri- 
marily because the IC limits are all roughly comparable 
as a result of the uniformity of the BAT sensitivity. If 
we instead take the best-fit IC normalizations, the non- 
thermal pressure of nearly all of the clusters is less than 
10% of thermal pressure with massive clusters typically 
at less than a few percent. Except for the least massive 
clusters in HIFLUGCS, the comparable magnetic field 
pressure assuming B = 1/iG corresponds to ~ 1% of the 
thermal pressure. Although we cannot measure B fields 
this large in the radio clusters, as long as B < 3-7/zG, 
it seems unlikely that non-thermal pressure support will 
significantly contribute in cluster centers, especially for 
the more massive clusters that are typically used for cos- 
mological parameter estimation. 



6.4. Future Outlook 

Can the survey observations with the BAT be im- 
proved, beyond the increase in sensitivity which comes 
with longer accumulating exposures? Perhaps the clear- 
est way forward is to better calibrate the spectral re- 
sponse of the BAT in narrower channel so that the fluxes 
are more reliable for steep thermal emission in the 14-24 
keV energy range. At present, we may be underestimat- 
ing source fluxes in these bands. If the first band is low 
by ~ 2ct and the second by ~ la, as suggested by the 
residuals in Figure 1 1 H our non-thermal limits will in- 
crease by about la — a small but non-negligible amount. 
The most straightforward fix is to remake the survey us- 
ing the BAT's native 80 channels instead of binning them 
into 8 channels that arc broad enough to be biased by 
the flux calibration with the Crab. With such improved 
data, this study can be repeated with a sample of all the 
known radio halo and relic clusters to definitively detect 
the non-thermal excess hinted at in the stacked "Radio" 
subsample considered here, if it exists. 

Ultimately, any IC detections, especially if marginal, 
will have to be confirmed by the upcoming missions with 
focussing hard X-ray telescopes, namely NuSTAI^ and 
Astro- -ffj- By resolving both contaminating point sources 
and the location of the hottest gas, these missions have 
the potential to achieve higher sensitivities than have 
thus far been possible. Also, these telescopes' narrow 
FOVs are well-suited to the smaller angular extents of 
mini-halos, making them ideal targets even though the 
radio flux densities are typically lower. 

We particularly owe the Swift BAT team a hearty 
thanks for uniformly processing the tens of thousands 
of individual pointings that make up the BAT survey 
that allow it to be such an incredibly useful resource 
for studies like ours. DRW and CLS were supported in 
part by NASA through Suzaku grants NNX08AZ99G, 
NNX09AH25G, and NNX09AH74G, and XMM-Newton 
grants NNX08AZ34G and NNX08AW83G. YYZ ac- 
knowledges support from the German BMBF through 
the Verbundforschung under grant No. 50 OR 1005. 
The XMM-Newton project is an ESA Science Mission 
with instruments and contributions directly funded by 
ESA Member States and the USA (NASA), and it is 
supported by the Bundesministerium fur Wirtschaft und 
Technologie/Deutsches Zentrum fur Luft- und Raum- 
fahrt (BMWI/DLR, FKZ 50 OX 0001) and the Max- 
Planck Society. Basic research in radio astronomy at 
the Naval Research Laboratory is supported by 6.1 Base 
funding. This research was also supported in part by 
an appointment to the NASA Postdoctoral Program at 
the Goddard Space Flight Center, administered by Oak 
Ridge Associated Universities through a contract with 
NASA. 
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Component 1 




Component 2 
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a Parameters for the 2T model are unconstrained. 

b Normalization of the APEC thermal spectrum, which is given by {lCP 14 /[47r(l + z) 2 D 2 A \} J n e riH dV, where z is the redshift, 
Da is the angular diameter distance, n e is the electron density, uh is the ionized hydrogen density, and V is the volume of the 
cluster. 

c Value is the normalization of the power-law component for the T+IC model, which is the photon flux at a photon energy 
of 1 keV in units of photons cm -2 s -1 keV -1 . For the 2T model, the value is the normalization of the second APEC thermal 
model in units of cm -5 . 

TABLE 6 

Joint Thermal and non-thermal Fits to the EPIC and BAT Spectra 

(3-195 keV) 
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2T 
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0.299 
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i O1+0.36 
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IT 
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2T 


2.80 
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0.0150 


3.62 
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IT 
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564.78/730 
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ft Q C 
D.OO 
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o.iyo 


U.Ulol 


8.95 
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T+IC 


8 41 + - 44 
°- 4i -0.36 
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IT 
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631.47/766 




2T 


4.35 


0.388 


0.0519 


5.03 
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T+IC 


4 5Q+O.I6 
^• oy -0.12 


U.OOD_q 021 


0806+ 0013 
U.UOUD_ 0023 
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UOU .lO/t uu 


A3376 


IT 


|_ 1 19 
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0.0086l°°°J° 






63.94/84 




2T 


1.98 


0.523 


0.0025 


6.83 


0.0070 


1^ Q Q 1 / OO 

Oo.oT/oz 




T-UTP 
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°- zo -1.31 
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IT 
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U.UUoo 


6.87 


0.0112 
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T+IC 
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D - ' -0.56 
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U.OOO_ o 8 3 
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A3395s 


IT 


„ li oq 
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56.82/124 
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4.ol 


U.zoo 


U.UUoy 


7.64 


0.0039 
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T+IC 
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u - ZDy -0.134 
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R1504 


IT 
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940.16/788 




2T 


6.93 


0.412 


0.0198 


9.48 


0.0299 


OQO QQ iTQd 

yoy.007 /OO 




T-UTP 

_L Tib 
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n 444+ 094 


n47^+ - 0014 
u ' w '°-0 .0043 
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A0576 


IT 
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„ if) <\Uti 

u.ooo_ 083 


^ „ „ ^ -Mi nn9f; 
0.02311°°^ 






91.85/129 




2T 


2.72 


0.430 


0.0136 


6.37 


0.0113 


on a. k / 1 07 
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T+IC 


4 16+°' 74 


402+ 125 
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0216+ 0037 
U.UZ10_ 0039 
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A0754 


IT 


9 - 46 -0.42 
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0.2So_ 033 


n n«n«+0.0011 
0.0696_ 00 ii 






528.60/645 




2T 


8.13 


0.288 


0.0564 


16.61 


0.0143 


529.13/643 




T+IC 


9 40 +0.43 


n 904+0. 038 
U.^84_ n 033 


n ocnc+0.0011 
0.0695I 00 43 




< 0.0017 


531.98/644 
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IT 




n 9s9 +0.026 
U - /OZ -0.026 


0414+ 0014 
U.U414_ 0013 






330.74/443 




2T 


3.89 


0.297 


0.0343 


6.84 


0.0079 


330.58/441 




T+IC 


4 on+0.25 
-0.24 


29 o+0.038 
u.zao_ 035 


0399+ 0028 
u.uo»3_ 0031 




< 0.0016 


330.55/442 


A1060 


IT 


40+0.09 


n 9 S -?+ ' 025 
U.OOO_ 023 


0558+ 0016 
u.uooo_ 0012 






519.24/641 




2T 


3.43 


0.383 


0.0290 


3.42 


0.0268 


519.24/639 



Wik et al. 



TABLE 6 — Continued 









Component 1 




Component 2 








kT 


abund 


Norm. 


kT 


Norin. c 
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Model a 


(kcV) 


z© 


(cm- 5 ) 


(keV) 




X 2 /dof 




T+IC 
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IT 
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2T 
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T4.Tr 1 
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2T 
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0.0203 


Z ( o.lo/o ( U 
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0.0019l°;°°?| 
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2T 
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IT 
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°- '°-0.04 


„ un m 3 
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2T 
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3.78 
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A1644 


IT 
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0.0205 


6.85 
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„ „ , _i_n (-,') 

6-94±g;ll 


„ , , _ _i_n Ofis 
u - 4il -0.064 
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u.uoo±_ 0013 
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2T 
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IT 
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2T 
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0.249 
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2T 


1.43 


0.688 
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801 + 1 701 
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IT 
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2T 
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2T 


5.39 
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0.0081 
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0.0087 
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T+IC 
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IT 
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o.iosgt 9 ;^ 1 ! 3 
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2T 
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/.oy 
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T+IC 
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IT 
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0.07781°;°°°! 
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2T 
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< 0.0009 


1 98 62 /28fi 
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IT 
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2T 
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n ^84+°- 321 
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2T 
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IT 
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420.99/524 




2T 
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IT 
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2T 
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IT 
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Component 1 Component 2 







kT 


abund 


Norm. 


kT 
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V KLV ) 


i© 
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2T 


4.42 


0.344 


0.0191 


4.65 


0.0168 


429.73/517 




T+IC 
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429.76/518 
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IT 
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2T 
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U.UUdl_ 0021 
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IT 
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2T 


5.55 


0.238 


0.0234 
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0.0147 


56.43/97 




T+IC 
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n nqR^+ 0041 

U.UOOO„g 0043 
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IT 


4 5O+0-14 
^• oy -0.13 
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573.33/746 




2T 


4.39 


0.367 


0.0498 


4.80 


0.0496 
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T+IC 
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< 0.0017 
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IT 
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2T 
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366.89/505 
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0442+ 0021 
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IT 
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2T 
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0.330 


0.0213 


7.25 


0.0286 
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T+IC 


7 70+ ' 47 
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0498+ 0014 
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IT 
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0242+ ' 0017 
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2T 
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0.0134 
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55.55/117 




T+IC 
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IT 
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2T 
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1156.48/1319 
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IT 
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2T 
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0.0166 


3.90 


0.0088 


212.44/281 




T+IC 
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IT 
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2T 


3.41 


0.622 


0.0186 
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0.0016 
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IT 
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u.uzoo_ 0009 






310.13/407 




2T 


3.43 


0.317 


0.0124 


4.39 


0.0114 


310.12/405 




T+IC 


o 7s +0.29 
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n n997+ 0020 
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< 0.0012 
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IT 
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62.69/79 




2T 


0.44 


0.297 
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5.76 


0.0153 


61.04/77 




T+IC 




0.269±°;i« 


0179+ 0030 




< 0.0008 


62.73/78 


A2657 


IT 


c 7q+0.64 
o - ,,:> -0.60 


rion +0.068 
U.ZO±_g Q66 


n242+ 0016 






176.95/235 




2T 


2.91 


0.267 


0.0054 


6.45 


0.0199 


177.03/233 




T+IC 


t- --+0.60 
°- ' ' -0.64 


n 9 c 9 +0.067 
U -^ J ^-0.068 


0241 + 0017 

U.UZ41_ 0Q21 




< 0.0007 


176.95/234 


A4038 


IT 


, og + 0.12 


Q45+O.O26 
U.J*O_ 02S 


0559+ 0019 






579.52/726 




2T 


3.07 


0.350 


0.0298 


3.73 


0.0267 


579.38/724 




T+IC 


3 36+ 014 
• 3 -"-' D -0.18 


n qt-9+0.045 
U.OOZ_q 031 


0554+ 0024 
u.uoo4_ 0032 




< 0.0013 


579.42/725 


A4059 


IT 


4 4Q+0.23 


0.428l° ;°| D 


0329+" 0012 
u.uj^y_ 0011 






259.39/427 




2T 


4.39 


0.431 


0.0162 


4.39 


0.0168 


261.26/425 




T+IC 


4 44+ .' 2 . 2 . 

^•^-0.22 


495+ - 040 
U -^ ZO -0.033 


0329+ 0013 
u.uc>zy_ 0016 




< 0.0005 


259.32/426 



a Parameters for the 2T model are unconstrained. 

b Normalization of the APEC thermal spectrum, which is given by {10 -14 /[47r(l + z) 2 D A ]} f n e nn dV, where z is the redshift, 
Da is the angular diameter distance, n e is the electron density, uh is the ionized hydrogen density, and V is the volume of the 
cluster. 

c Value is the normalization of the power-law component for the T+IC model, which is the photon flux at a photon energy 
of 1 keV in units of photons cm~ 2 s _1 keV -1 . For the 2T model, the value is the normalization of the second APEC thermal 
model in units of cm -5 . 



